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ABSTRACT 



A mean streamline analysis of the effect of stator blade orienta- 
tion on the performance of an axial flow compressor was performed by 
means of a computer program. Measurements were made on a 3-stage axial 
flow compressor at the Naval Postgraduate School at six stator stagger 
angles between 23.8^ and 44.3^ for a fixed orientation of the rotor 
blades. Maximum efficiency and pressure ratio were measured at a 
stator stagger angle of 31.8^. Results at other blade settings showed 
that by varying stator stagger angle with flow rate optimum efficien- 
cies and pressure ratios can be achieved over a wide range of operating 
conditions . 

The results of the analysis were compared with the measured results. 
Suggestions are made for improving the manner of adapting cascade test 
data to performance predictions. 

By applying a non-dimensional deflection coefficient it could be 
shown that minimum work input corresponded to maximum efficiency. 

The test compressor has a tip diameter of 36 in. and a hub/tip 
ratio of 0.6. The blading tested is of the free-vortex type with a 
design degree of reaction of 0.5. Tip speed was about 185 ft/sec. 
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SECTION 1 



INTRODUCTION 



Axial flow compressors have a wide range of application. They 
are usually matched to a turbine in a set. Both open-cycle sets, 
such as aircraft turbo-jet engines, and closed-cycle sets for power 
generation are in service. In every application compressor efficiency 
is the critical factor in set performance. Accurate predictions of 
off -design compressor performance must be made during the design stage 
to be able to make corrections prior to manufacturing. 

Theoretical analysis of the flow in an axial compressor involves 
three-dimensional partial differential equations of a complex nature (1). 
In these equations it is difficult to account for real gas effects such 
as boundary layer growth on machine walls and blade surfaces c The 
resulting wakes behind blade rows create non-uniform conditions; hence, 
it becomes necessary to base prediction methods on experimental data. 

These data would be obtained best on an actual compressor. How- 
ever, most compressors have relatively short blades and small flow 
annuli. Even if only small pressure probes were inserted between the 
rows of blades, the flow in these machines would be substantially 
altered. Test rigs must therefore be used that have large dimensions. 

A first approximation of the flow in axial turbomachines can be 
obtained in rectilinear cascade test rigs. The intersections of a 
co-axial cylinder with the blades of a row produce a series of iden- 
tical and identically oriented profiles which are unwrapped into a 
plane to establish the corresponding rectilinear cascade. A finite 
number of straight blades with profiles similar to, but larger than 
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those obtained, is then arranged in a rectangular duct through which 
air is blown at the appropriate inflow angle. Flow surveys taken ahead 
of and behind the blade row give experimental data of air turning angle 
and total pressure loss for various blade geometries. Such a recti- 
linear cascade is in use at the Naval Postgraduate School. The results 
of extensive testing of compressor airfoil shapes at NASA have been 
summarized by Lieblein (2). 

The usual problem in compressor design consists in selecting flow 
areas, blade shapes and blade layout to satisfy design specifications. 
Additional studies are frequently necessary if the design performance 
is not reached. An example is the study of Vavra which was made when 
the so-called Clark CSN-1 compressor of the ML-1 nuclear gas turbine 
of Aerojet-General Nucleonics failed to perform as required. Vavra 
analyzed the design changes proposed by the manufacturers and made 
further recommendations (3). However the ML-1 project was cancelled 
before the improved compressor could be built so that it was not 
possible to verify the suggested changes by experiments. 

The so-called inverse problem consists in predicting the perform- 
ance of an existing machine with a known blading by means of rectilinear 
cascade data and other experience factors. Such an approach was carried 
out by Gibbons and Bartels (4) for the 12-stage Allis-Chalmers axial 
flow compressor which supplies air to the turbine test facilities at 
the Naval Postgraduate School. Because of the unorthodox performance 
of the first stages, and also due to small flow annuli and short blade 
heights, it was not possible to predict the performance of the machine 
accurately. However, the analysis gave important indications for 
design changes to improve the performance of the compressor. 
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This thesis makes another attempt to solve the inverse problem 
for a 3-stage axial flow compressor built by the California Institute 
of Technology with the support of the Office of Naval Research. The 
compressor is sufficiently large so that the insertion of pressure 
probes causes relatively small flow perturbations. The machine was 
designed to permit variations in blade shapes, blade angles, tip 
clearances, and staging. The installation at the California Institute 
of Technology is described by Bowen, ^ (5) This report also 

describes a theory of perfect fluid flow in axial flow turbomachines. 
The results of this theory are compared with the measured data of the 
first stage of a so-called **free-vortex** type blading. Part 2 of the 
report by Bowen, (b) , gave a detailed investigation of multi- 

stage flow for both **free-vortex** blading and a more complex ** solid 
body rotation” blading. Measurements of blade skin friction losses 
were greater than expected from cascade tests. The growth of boundary 
layer along the inner and the outer annulus walls of the flow was less 
than expectedc Alsworth and lura (7) conducted extensive and detailed 
measurements of flow patterns in a single stage of ”free-vortex” 
blading. They carried out accurate measurements of the blade skin 
friction losses and the radial distribution of work input. The final 
report from the California Institute of Technology was a hot-wire 
anemometer study of compressor stall by lura and Rannie (8). They 
observed that stalled flow regions rotate in the direction of blade 
rotation without changing shape but with a speed that is not propor- 
tional to rotor speed. 

After the Turbo-Propulsion Laboratory of the Naval Postgraduate 
School was built, the compressor was relocated there. Since then 
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it has been used for laboratory courses to supplement instruction 
in basic theories of turbomachines. 

The present study was conducted for this compressor because of the 
possibility of changing its blade angles. For fixed stator and rotor 
blade angles a highly peaked curve of efficiency versus flow rate is 
usually obtained with axial flow compressors. At a slightly changed 
stator blade angle the peak will be displaced somewhat. A family of 
efficiency curves for various stator stagger angles is expected to 
yield a flatter efficiency curve over a wider range of flow rates for 
fixed rotor blade angles. The wide variation in flow rates required 
of the engines for the supersonic transport has led to design pro- 
posals by General Electric for compressors where the stator blade 
angles can be changed during operation of the jet engine. Hence, the 
thesis topic is of current interest. 

For the tests it was necessary to obtain extremely accurate 
pressure measurements. An extremely sensitive bourdon tube pressure 
indicator was added to the instrumentation. A permanent Pitot-static 
tube was installed, and accurate flow rate calibrations were carried 
out. A parallel theoretical effort produced predictions of com- 
pressor performance by an existing computer program which was modified 
and adapted for use on the I. B. M. Model 360 computer system installed 
at the Naval Postgraduate School. 

The author wishes to express a debt of gratitude to the faculty 
of the Department of Aeronautics of the Naval Postgraduate School for 
his engineering education. Particular thanks are due to Professor 
M, H. Vavra for his enthusiasm and guidance in showing what can be 
done with that education. 
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SECTION 2 



0. N. R. 3-STAGE AXIAL FLOW COMPRESSOR 

The compressor installation is shown in Fig. 1. The compressor 
inlet is in the background of the picture. The exit section is 
equipped with a throttle valve. A torque meter is installed on the 
shaft between the drive motor and the compressor. The overall dimen- 
sions of the equipment are given in Fig. 2. 

A basic criterion of the design of the compressor was flexibility 
of operation. Each blade row may be removed. The angle setting of 
each blade is adjustable in 0.5^ increments. In order to minimize 
flow disturbances by pressure probes the outside diameter of the blad- 
ing is 36 in. For the same reason^ the compressor has a hub/tip ratio 
of 0.6, giving an inside diameter of 21.6 in. and a blade height of 
7.2 in. To permit simulation of conditions in a multistage unit, the 
machine has three stages. One row of inlet guide vanes simulates the 
effects of previous stages. Two rows of exit guide vanes remove the 
whirl component after the third stage. To insure that general align- 
ment and tip clearances are maintained the compressor casing is very 
rigid. 

The outer casing consists of two cast iron half-cylindrical 
shells bolted together along the horizontal plane through the axis. 

The rows of inlet guide vanes, stator blades, and exit guide vanes 
are held by bolts extending through the casing (Fig. 3). Inside the 
casing a hollow steel shaft carries three cast iron drums for each 
of the rows of rotor blades. The shaft is supported by two roller 
bearings whose outer races are pressed into the supporting strut 
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assemblies (Fig. 4). These assemblies rest in close fitting channels 
in the casing (Fig. 5). The six assembly struts have symmetrical air- 
foil sections. Figure 6 shows the nine blade rows after assembly. 

The casing was designed to give maximum accessibility for mea- 
suring instruments. Six rectangular instrument ports are located in 
the upper half at 30^ from either side of the vertical (Fig. 3). They 
are placed in the second and third rotor planes, the first, second and 
third stator planes, and behind the third stator. The ports hold a 
special instrument carriage permitting detailed flow surveys in any 
axial plane. The casing has numerous radial survey holes whose loca- 
tions are specified in Fig. 7. 

The blades are made of ALCOA 356 aluminum alloy without heat treat- 
ment. There are thirty rotor blades and thirty-two stator blades per 
row. The details of construction of the blades are illustrated in 
Figs. 8 and 9. The blading is designed to have a degree of reaction 
of 0.5 at the mean radius. The rotor blades are twisted by 49^ from 
hub to tip, and the stator blades by 13°. At the mean radius, the 
important blade parameters are: 



Parameter 


Rotor 


Stator 


Camber angle 


20.32° 


29.98° 


Design stagger angle 


43.80° 


28.80° 


Thickness- to-chord ratio 


0.1 


0.1 


Chord length 


2.60 in . 


2.60 in. 



The values of the blade parameters at other radii, and the method 
used for calculating the thickness distribution, are specified in 
Ref. 5. 
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The point of maximum thickness was set at 0.35 of the chord 
length from the leading edge. The thickness was modified over the 
rear 15 per cent of the section to provide a trailing edge thickness 
0.02 in. The thickness distribution is applied about a parabolic 
mean camber line. Figure 10 is a picture of a ‘'free-vortex** stator 
blade. Tip clearances of 0.020 in. for stator, and 0.037 in. for 
rotor, are maintained by the use of shims at the point of attachment. 

The attaching device permits variation of rotor and stator blade 
stagger angle. Figure 11 shows how these changes are made. The blade 
shaft extends through the rotor drum at the hub. An adjusting plate 
is secured to the blade shaft by a tapered pin. The plate has a series 
of eleven taps on a circle about the shaft axis which are spaced at 
intervals of 4.5^. Between the adjusting plate and the rotor drum is 
a fixed sector which is attached to the drum by a pin which guarantees 
its alignment. The fixed sector has eleven holes which are spaced at 
intervals of 4.0^ of arc about the blade shaft axis. Blade alignment 
is maintained by a set screw through both holes in the adjusting plate 
and the fixed sector. If the center holes of the plate and the sector 

are lined up, the blade is at the design stagger angle of 43.8^ for 

the ** free- vortex’* rotor. An angle change of 4.0^ can be accomplished 
by keeping the set screw in the center hole of the adjusting plate and 
moving the blade so the screw is inserted into the next hole of the 
fixed sector. Small angle changes of 0.5^ are accomplished by moving 
the blade slightly so that the two holes immediately adjacent to the 
center holes are lined up since blade angle is changed only by the 
difference in arc between the two holes. A similar arrangement is used 

for the stator blade settings (Fig. 10). Figure 12 illustrates stator 

blade settings which are one degree smaller than the design stagger 
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angle. The two blades on the right-hand side of the figure belong to 
the first stator row. The set screw is two holes away from the center 
or reference holes. The blades on the left-hand side of Fig. 12 belong 
to the row of inlet guide vanes. They are set at an angle which is by 
0.5^ smaller than the design stagger angle. 

The assessories to the basic compressor will be described by fol- 
lowing a path from the inlet to the exit and to the power source. The 
instrumentation of the compressor will be described in its appropriate 
place in the same sequence. Details of instrument calibrations are 
contained in Appendix A. 

The inlet duct is seen in the background of Fig. 1. It consists 
of a screen, an entrance bellmouth, and a length of straight pipe. 

Figure 13 shows the large mesh screen which prevents the ingestion of 
foreign matter from the compressor bay apron. Mounted on the screen 
is a mercury thermometer used to determine ambient temperature. It 
has provisions for psychroraetric analysis of the incoming airo 
Immediately behind the screen is the bellmouth. It changes the dia- 
meter rather abruptly from almost two diameters at the flare, to 
36 in. in the axial distance of about 15 in. A cylindrical inlet duct 
two diameters long connects the bellmouth to the compressor proper. 

The inner surface of the duct is enameled to give smooth flow surfaces. 

The flow rate through the compressor was determined by surveys in 
the inlet duct. Details of the calibration are given in Appendix B. 

The survey plane is about midway between the bellmouth and the com- 
pressor. The surveys were made with a Prandtl-type Pitot tube. In 
Fig. 14 the probe is installed for taking a vertical traverse. It is 
possible to disassemble the probe while the compressor is running (Fig. 15). 
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It may then be remounted for a horizontal traverse. Traverses at 30^ 
and 60° from the vertical direction are possible. These traverses can 
be made by unbolting the entire inlet duct from the compressor and 
rotating it on its cradle. The design of the probe prevents measure- 
ments closer than 0.25 in. to the inlet duct wall. 

A honeycomb straightener is installed behind the inlet Pitot tube 
to equalize the flow entering the compressor (Fig. 14). Attached to 
the forward strut assembly is an ogival wooden fairing (Fig. 5) to 
provide a smooth transition of the flow from the inlet duct to the 
compressor annulus. 

Pressure measurements may be taken at any one of the radial 
survey ports described. Figure 16 shows the probe used for these 
measurements. It is a United Sensor three-hole probe, model YC-120. 

A central hole measures stagnation or total pressure. A static port is 
located on each of two faces of a wedge at an angle of 45° from the 
total pressure hole. The static pressures are balanced on a water 
manometer to insure alignment of the probe in flow direction. The 
probe is mounted on a probe holder which has vernier scales permitting 
radius adjustments to 0.01 in. and angle adjustments to 0.1°. 

The accurate determination of flow rate is essential to compressor 
analysis. For this purpose a modified Prandtl type Pitot-static probe 
was installed. Figure 17 shows the installation details. The probe 
is inserted in a radial survey hole located ahead of the inlet guide 
vane row. It protrudes forward into a channel between two adjacent 
struts of the supporting assembly which is shown in Fig. 5. The probe 
is approximately at the center of the channel between struts. The 
probe is positioned in the radial survey hole by a brass plug, and 
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secured to the casing by a simple locking device. Figure 18 shows the 
probe and its locking device in the background. In the middle of the 
figure are seen the three-hole probe and its holder located between the 
inlet guide vanes and the first rotor. In the foreground, a plug, 
originally located at the same axial position as the probe, has been 
removed to show a hole through which the probe can be inserted. 

Pressures obtained by the probes were measured by a Texas Instru- 
ments Fused Quartz Precision Pressure Gauge which is shown in Fig. 19. 
The velocity head in the inlet is about 0o6 in. of water. A low- 
pressure bourdon tube was obtained for the gauge which measures from 
0 to 166 in. of water with the instrument reading from 0 to 200,000 
counts. The calibration constant of the low pressure bourdon tube was 
determined to be 240*423 counts/psf of the pressure gauge (See Appendix 
A.). The pressure gauge operates both in manual and servo modes, the 
latter providing automatic nulling of the unit. 

Since flexibility of pressure selection was desirable, two 
Giannini Sp-lOlA pressure scanners were connected to the pressure and 
reference ports of the bourdon tubeo One of the scanners is seen in 
the foreground of Fig. 19. The arrangement of pressures to the 12- 
channel switches is shown in Fig. 20. Later in the study a simpler 
system of five valves and two manifolds was constructed (Fig. 21) to 
eliminate leakage flows that seem to have occurred in the scanners. 

At this time the original 3-hole probe was placed permanently behind 
the third stator (Fig. 22). Another 3-hole probe was placed ahead of 
the first stage rotor. This arrangement permitted direct measurement 
of the pressure increase across the three stages^Z^ P^ g. 
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Figure 22 also shows the instrument traverse carriage mentioned 
before. It was used for a survey of total pressure in peripheral 
direction behind the third stator at two different radii (See Appendix 
C.)* it is possible to vary the radial locations of the probe with an 
accuracy of 0.01 ino The probe may be rotated through 360° about its 
axis. The carriage can move the probe 15° in peripheral direction^ 
which covers a whole blade spacing since a stator blade channel covers 
11.25°^ and a rotor blade channel 12.0° in peripheral direction. The 
locations of the radial survey holes with reference to the indicated 
meridional angle of the carriage are shown in Fig. 23. The represented 
stator profiles are at the design stagger angle at the mean radius. 

The actual stagger angle measured with respect to a plane through the 
axis is the complement of the indicated angle. 

A short cylindrical duct connects the compressor to the exit 
elbow, which is equipped with sheet metal turning vanes to minimize 
losses. A transition piece connects the elbow to the throttle valve 
which is seen in Fig. 1. The valve consists of two rectangular metal 
doors which move in slides. The position of the doors is controlled by 
a right- and left-hand threaded lead screw which is rotated by a 
Graham variable-speed transmission. A revolution counter indicates 
the approximate valve opening. 

The compressor is driven by a 50 HP Fairbanks Morse induction 
motor. It requires a three-phase, sixty-cycle , 440 volt power supply, 
and can operate at two fixed speeds; namely, at about 900 rpm or 
1200 rpm. 

Attached to the motor shaft is a Baldwin-Lima -Hamilton SR-4 
torque meter type A. It is shown in Fig. 24 during calibration with 
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known weights attached to a lever. In the background may be seen a 
Brown Instruments strain gauge readout which is of the standard 
Wheatstone bridge type. A conversion constant of 4.1565 ft-lb/lb of 
torque meter readout was determined (See Appendix A.). Under the 
protective cover next to the torque meter in Fig. 24 is installed a 
six lobe flux cutter to measure the motor? speed by means of an 
electronic counter. 
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SECTION 3 



FLOW RATE CALIBRATION 



The permanent Pitot-static tube was installed to provide rapid 
and accurate measurements of volume flow rate. At a given throttle 
valve setting this probe was calibrated against the result of two 
traverses of the cylindrical inlet duct. Through-flow velocities in 
the flow annulus ahead of the inlet guide vanes were varied from 
110 ft/sec to 70 ft/sec, so that compressibility effects could be 
ignored. The objective of the calibration was to determine a cali- 
bration constant which gave the volume flow rate when multiplied by 
the velocity measured at the permanent probe. 

An initial survey of the inlet duct was made at the highest pos- 

3 

sible flow rate of 512 ft /sec. Figure 25 shows a plot of the mea- 
sured velocities of both horizontal and vertical traverses. The 
smooth velocity profile in the duct is evident. Figure 26 is a graph 

of the measured velocities in the inlet duct at the lower flow rate of 
3 

310 ft /sec, which shows a higher degree of scatter. For these tests 
the pressure indicator was operated in the servo mode. Even in the 
meter mode the indicated pressure oscillated considerably. Pressure 
fluctuations seemed to have been initiated by the vibration of the 
brass rod on which the traversing Pi tot- static tube was mounted in the 
inlet duct. The pressure indicator corrects a disparity between the 
angular position of a mirror of the bourdon tube and the indicator 
dial at the rate of the full scale reading of 200,000 counts in 120 
seconds. With the oscillating pressure applied to the sensitive 
instrument, oscillations of 50 counts were observed for a pressure 
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corresponding to a reading of 600 counts. To reduce these oscillations, 
variable length capillary tube damping devices were installed both in 
the total and static pressure lines of the traverse Pitot-static tube. 



Details of the arrangement and the results of tests are shown in 
Appendix B-1. 

A special survey was made to determine the thickness of the bound- 
ary layer on the walls of the cylindrical inlet duct. However, con- 
struction of the traverse probe prevented taking readings closer than 
0.25 in. from the wall. Figure 27 is a plot of the readings at high 
and low flow rates. In each case it is evident that large losses occur 
near the walls. With the velocity distributions of Fig. 27 the flow 
rate is about 2.5 per cent smaller than the value obtained without 
considering the changes near the walls. 

Data were taken in the boundary layer at intervals of 0.25 in. 
from 17.75 in. to 16.00 in. In the main stream^ data were taken at 
4.0 in. intervals from 16.00 in. to the centerline. With horizontal 
and vertical traverses, four points were obtained at each radius and 
two at the axis for a total of 46 locations. At each station the 
following measurements were taken: 

- velocity head in the inlet duct (counts) 

- velocity head of the permanent Pitot-static probe (counts) 

- static pressure with reference to atmosphere in inlet 
duct (counts) 

- static pressure with reference to atmosphere at permanent 
Pitot-static probe (counts) 

- temperature at the inlet screen (^F) 

- barometric pressure (in. of Hg) 






si 



so 



Bar 



Bar 



- temperature of column of mercury ( F) 
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The barometric pressure P (in. Hg) was converted into pounds 

Bar 



per square foot to obtain the atmospheric pressure P by 

A 

(psfa) 



( 1 ) 



Pf. = Pear WI. 4 G 7 ) 

where the constant was obtained by correcting the height of the column 
of mercury for temperature variations in specific gravity of the liquid. 
The inlet temperature was converted to absolute temperature, and was 
considered to be the total temperature T^^ because of the low veloc- 



ities, or 



(°R) 



( 2 ) 



The static pressure at the inlet was converted from the meter 



(psfa) 



(3) 



reading to an absolute pressure by 

P' = p — — 

^ ^ ^40.4^3 

where the constant was obtained by the calibration procedure described 
in Appendix A. The pressure P^ measured at the permanent probe was 
obtained in a similar manner. The local air density p, at the inlet 






was computed from 



^ O.OOZlya Pj 51^ ^ (Ib sec^/ftS 



(4) 



Cl4.C.%)(l44) \ 



T 



•to 



and similarly at the permanent probe. The velocity head in the 
inlet q^^ was converted from the meter reading. 



%ii = 



240. 423 



(psf) 



(5) 



The dynamic head at the permanent probe q^^ was obtained similarly. 



The velocity in the inlet is then 



V 



0 w 









(ft/sec) 



( 6 ) 
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A similar procedure gave the velocity at the permanent probe. The 
arithmetic average of the 46 values of obtained during the carrying 
out of the traverses gives the average value of 



4 G 



V 



od V 



= — 7 V • 

46 V 



(ft/sec) 



J^l 



The measured values of were adjusted to account for variations in 



flow rate due to slight compressor speed variations 



v,.,d = V,. 



V, 



(ft/sec) 



The four values of obtained at each radius were averaged by 



Vr 



I a V 






(ft/sec) 






The volume flow rate was first calculated by neglecting the losses in 
the boundary layer to give <> In general the volume flow rate 

is obtained by 



V" 




(f t^/sec) 



If the above equation were integrated directly, the velocity at the 
axis would not be included in the volume flow rate calculation 
because the zero radius is a member of the product under the integral. 
Therefore the integral was replaced by a summation. 



V^t ' ^ (£tW) 

The quantity A A = 2 Tf r h/144 represented that part of the annulus 

area in the center of which V. was measured, where h is the height 

iav 

of the area. This method assumed that the velocity measured at the 
radius of 16.0 in. exists also at the wall. A volume flow rate 



36 



deficiency was then calculated to account for the difference in 
velocities at 16.0 in. radius and those measured in the boundary 



layer, or 




By a trapezoidal integration the volume flow deficiency becomes 







~z7f Vq r Jy' (ft^/sec) 



i(,.o 



Then the actual volume flow rate is 







A calibration constant CCC for each run was obtained from 



CCC - — 



3 

(ft /sec per ft/sec) 



oav 



A first series of calibration runs was invalidated when a leakage 
hole was discovered in the traversing Pitot-static tube. After repairs, 
the calibration constants CCC were obtained with a maximum relative 
error of 0.1 per cent. The final calibration constant was 4.4050 



nently installed Pitot-static probe. 

A data reduction computer program, ONRFLO, for these calibration 
procedures^was generated for the I. B. M. Model 360 computer. A 
listing of the program is included as Appendix B-2. Sample output 
data from the final run is given in Appendix B-3. 



3 

ft /sec per ft/sec velocity obtained from the readings of the perma- 
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SECTION 4 



MEASUREMENT OF COMPRESSOR PERFORMANCE 



The objective of the study was to determine the effect of changes 
of the stator blade orientation on compressor performanceo It was 
decided to restrict performance investigation to the three stages only; 
namely, from a station ahead of the first rotor to a station after the 
last stator, or between locations SP-2 and SP-8 of Figo 7. The effects 
of the inlet, the inlet guide vanes, the exit guide vanes, and the 
diffuser were not considered. 

For tests at a particular stator stagger angle the barometric 

pressure P was determined first. For various settings of the 
Bar 

throttle valve the following measurements were taken: 

T^ - ambient temperature (°F) 

P^^ - static pressure with reference to atmosphere at the 

permanent Pitot-static probe (counts) 

- velocity head at the permanent Pitot-static probe (counts) 
total pressure with reference to atmospheric between the 
inlet guide vanes and the first rotor (counts) 



tig 



A P^ g - total pressure difference across the three stages (counts) 
T^ - reading of torque meter (lb) 

N - compressor speed measured by electronic counter (rpm) 

With the relations listed on p. 35 these measuring data establish 
the following quantities: 

(psfa) (1) 

( 2 ) 









1 4 - 45 ‘‘ 5,7 

• i 



(°R) 





240.^2 3 



( 3 ) 



38 



(psfa) 



0.00^376 Po 5f8.7 
(14.CU)(|40 



(5,S2xio“) (lb sec^/ft^) 
Tto 



(^) 



% 



ol 



%o 

240.^Z3 




^ 4oi 



(psf) (5) 

(ft/sec) (6) 



The absolute total pressure ahead of the first stage, or after the 



inlet guide vane, P , , is 

tiga 

p == p - 

^ 240.423 



(psf) 



(7) 



The pressure difference across the three stages A P^-SA obtained 
from 



aP 



^ ^3-S 



1 



240. 4Z3 

i meter reading 

- C4.I5-(,5) 



(psfa) 



From the torque meter reading T^ the actual torque T^^ is 



4, 



(ft-lbf) 



where the constant was obtained by the calibration procedure described 
in Appendix A. The angular velocity is with the measured speed, N 



Ca> ^ 



zm N 

CO 



(rad/sec) 



Both volume and weight flow rates were calculated. The volume 



flow rate 






is then 



^ V, ( 4.4 050) (ft^/sec) (8) 

where the constant was obtained by the calibration procedure described 
in the preceding section and Appendix B. The weight flow rate w is 
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. _yp 32.H4 Pa 518,7 _ A. / 
w - vc , . . -p (.7. 88i>) (Ibm/sec) 



O+.t'ICX'^^) X 



( 9 ) 



to 



For comparison with the prediction program a referred flow rate w^ is 



defined by 



\Nyf. = 



W 






P. /l44 



(Ibml'^^R/sec psia) (10) 



Two basic indices of performance were computed; namely, the 
3-stage pressure ratio ^ , where 

o ^^-5A 



r 3-S 



P-fetqd 



( 11 ) 



and the 3-stage total- to-total efficiency ^ obtained from 



13-: 



^Wfh 



PwA 



(12) 



where the theoretical and actual horsepower, respec- 
tively. The theoretical power is given by 

AvT+h = W ,-5^778.3) (ft-lbf/sec) 



where is the isentropic specific enthalpy difference for a 

frictionless compression from . to P^, after the third stator« 
^ tiga td 



Assuming a perfect gas there is 
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/ 

3-5 



(Btu/lbm, R) 



where g i^ ^be isentropic temperature difference from the total 

temperature ahead of the first stage, to after the third 

stator (See Fig. 28.). Since very nearly equals the ambient 



temperature there is 
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Tt, 



IT'' ^ 

1 

T 

\ ^ '9 



(Btu/lbm, R) 



With the isentropic relation 



r. 



/ 



td 



-td 






T 



pf, 



9' 



the final expression for the determination of the theoretical power is 



fw+K = ^ (ft-lbf/sec) (13) 



The actual power is obtained from 

(ft-lbf/sec) 

The total-to- total efficiency at the three stages is then 



Pw/^ ” LO 



a 






3-S ■ T 



% 



OH 



[(?,«)■*■- 1] 



778.3 



(14) 



For comparison with other compressor data certain non-dimensional 

parameters are introduced. The average flow function is the 

ratio of through flow velocity V and peripheral speed at the mean 

a 

diameter U , or with 
av 

"" 451.4/144 "" Co.221) (ft/sec) 

and, for the mean diameter of 14.4 in. 




00 



14,4 

I Z 



there is 



(|2)(I44) 

CO 0+.4)C45|.^) 



(ft/sec) 



OO 



(0.18AZ') 



(15) 
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An average dimensionless stage pressure coefficient is defined by 

0 U ^ 

^ ^ Xc 



_ -q/ 

*■* * n 



^SJ 



CW 



(psf) 



or 



Y 



aP 



3-SA 



(16) 



3 U., 

The so-called dimensionless deflection coefficient is defined by 



aP. 
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(psf) 



Hence , 



X 



Y 



c\v 






^3-S 



(17) 



The performance calculations were made for the conditions along 
the mean radius of the stage^l4.4 in., by assuming that these conditions 
are representative of the performance of the 3-stage compressor. 
Detailed investigations which support this simplification are given in 
Appendix C. The data reduction was accomplished by the computer 
program ONRETAj which is listed also in Appendix C. 

Data were taken on 11 occasions for the following settings of 
stator stagger angle: 

Run 



1, 2, 3 
4, 5 
6, 7 
8, 9 
10 
11 



Stator Stagger Angle 
27,8° 

23.8° 

31.8° 

35.8° 

39.8° 

44.3° 
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SECTION 5 



PERFORMANCE PREDICTION PROGRAM 

The program used to predict the performance of multi-stage com- 
pressors is called AXC03. It is listed in Appendix D and represents 
an adaptation of a program by Vavra (3), modified by Gibbons and 
Bartels (4) . 

The method assumes that the flow is axisymmetric. Computations 
are made for the streamlines at the mean radius of a stage. A stage 
by stage analysis of the machine is performed. The output of one 
stage is used as input for the next stage. The computation is started 
at a flow rate less than that anticipated for surge, and the flow rate 
is then increased by small increments to cover the whole operating 
range of the machine. 

Program AXC03 uses two principal subroutines; namely, one called 
ROTOR and the other called STATOR. They establish the velocity 
triangles of a rotor or a stator row of blades, respectively, from the 
known geometry of the bladings, flow rate, and discharge conditions of 
the preceding row of blades. Both subroutines make use of three 
additional subroutines . 

Subroutine THEORY calcualates the flow conditions for minimum 
profile losses. Cascade data by Lieblein are used for this purpose (2). 
These data are presented as curves for zero-camber incidence angle 
^ ^ ^ ( O ’ slope factor n, zero camber deviation angle (€o) , and 

slope factor m, all given as functions of the inlet air angles 
or ^ j , with the blading solidity as a parameter. The solidity 
or is defined as the ratio of blade chord c and blade spacing s. 
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The curves described above have been expressed by polynomials in powers 
of solidity and inlet air angle which are used in subroutine THEORY. 

The accuracy of the polynomials has been verified, and the errors are 
listed in Appendix D. The incidence angle for minimum loss is deter- 
mined by an iteration procedure which changes the inlet air angle until 
agreement is reached within a tolerance of 0.02°. The deflection which 
occurs at the minimum profile loss can then be calculated with the 
established polynomials also. 

The subordinate subroutine CASCAD calculates the performance of 
the blading at the actual incidence or inlet flow angle of the blade 
row by using a graph of Lieblein that establishes the change of the 
deviation angle with incidence angle ( d S /d C ) as a function of 
inlet flow angle and solidity. 

In the earlier computer programs the curves of this graph were 
expressed by a polynomial that covered the range of solidities from 
0.0 to 1.8. Since the curves have nearly exponential character, which 
cannot be expressed with ease by polynomial functions, it was found 
that an error of the order of 10 per cent was possible. Therefore, 
a new polynomial was established for solidities between 0.6 and 1.4 
which limits the error to less than 1.9 per cent. 

Subroutine CASCAD further utilizes experimental limits on blading 
performance by Howell (9) which are illustrated in Fig. 29. Howell 
found that for most cascades the ratio of the profile drag coefficient 
Cj^ at an arbitrary incidence angle L and the minimum loss 

incidence is a unique function of the quantity > 

where 6 is the optimum flow deflection for the incidence angle 
Since the applicability of the data presented in Ref. 2 is limited to 
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design point calculations, where the incidence angles are close to the 
nominal incidence | q > the curves cannot be used for off-design analy- 
ses without restrictions. Without such limits the flow deflections 6 
could be increased simply by increasing the incidence angle i y since 
the data of Ref. 2 do not establish criteria for flow separations 
which are associated with a radical increase in profile losses. Exper- 
imental results of Howell show that the flow deflection 6 at incidence 
t should not exceed the optimum deflection € by more than 25 per cent 
to avoid excessive losses. He states that this limit occurs if the 
profile loss coefficient Cj^ is about twice the minimum loss coefficient 
^Dmin' solid curve in Fig. 29, labeled 6/6 is supposed to be 

the relation between ^ ^ /o obtained from subroutine 

CASCAD with the data from Ref. 2. It will be assumed that this curve 
holds for negative values of ( /o )/^ / ^ only. For positive values 
of this parameter the dashed curve will be applied to satisfy the 
limitations of Howell. This curve is somewhat arbitrary but represents 
the best estimate for the actual operating condition of compressor 
cascades. Figure 29 further shows a curve relating the ratio C^/C^ . 
to 9 which is adapted from Ref. 9. Both this curve and 

the dashed curve for the deflection ratio have been expressed analy- 
tically in subroutine CASCAD. If the ratio C^/C^ . exceeds a value of 
2.0 for positive quantities (L“L|q )/^|o > program prints out 

"surge in blade row xx" and proceeds with calculations at an increased 
flow rate. The positive value of at which this occurs is 

the input quantity SMAX. If exceeds a value of 2.0 at 

incidence angles L smaller than | o > the program prints out "minimum 
pressure in blade row xx" and halts computation. The negative value of 
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at which this occurs is the input quantity SMI. In 
both instances the symbol xx refers to the number of the blade row, 
and indications are given also whether the row is a stator or a rotor. 
For any incidence angle within the useful operating range, the values 
of the deflection ratio and of the ratio C^/C^ . are computed by sub- 
routine CASCAD and control returned to the blade row subroutines 
ROTOR or STATOR. 

The efficiency of a blade row is determined by subroutine ETACAL. 
Losses due to tip clearances, secondary flows, and wall friction are 
taken into account with the relations proposed by Vavra (l)o 

If conditions of surge exist, the flow rate is increased by about 
1.0 per cent and the entire computation is repeated. If the incidence 
angles are in the so-called useful operating region of Fig. 29, the 
output of subroutine ROTOR or STATOR becomes the input for the sub- 
routine that calculates the next row of blades. Interstage data are 
printed out as the program progresses from one stage to the next. If 
all blade rows of a machine have been processed, the overall compressor 
efficiency and pressure ratio are computed for the particular flow rate. 
This flow rate is then increased by a specified amount, and the calcu- 
lating process is repeated. If the so-called minimum pressure is 
reached, the computation stops and a summary of the overall performance 
parameters is printed out. 

The details of the calculation of the flow through the inlet duct, 
the inlet guide vanes, the exit guide vanes, the diffuser and the dis- 
charge passages are not described in this thesis since it is concerned 
only with the performance prediction of the three stages of the com- 
pressor from a station ahead of the first rotor to the discharge at 
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the third stator. The methods applied for the flow analysis in these 
passages may be found in Vavra (3) and Gibbons and Bartels (4). 

According to Howell there occur additional effects that influence 
the off -design calculations. The growth of boundary layers on the 
walls of the annular flow channel will change the velocity profiles 
and can reduce significantly the effective flow area. The program 
takes account of this effect by a blockage factor that can vary from 
stage to stage. This boundary layer growth is responsible also for 
increasing peaks in the radial distribution of the axial velocity 
components in succesive stages. These higher velocities outside of 
the wall boundary layers produce smaller incidence angles for the 
main portion of the flow, thereby decreasing the actual flow deflection 
and reducing the work absorbed by the fluid in the stage. Hence 
in a compressor that consists of stages with identical bladings, 
the last stages will produce smaller pressure ratios than the stages 
at the compressor inlet. Experience shows the efficiency in successive 
stages is not reduced by the peaking of the velocity profiles, and 
their effects are usually taken into account by a so-called work-done 
factor that is about unity for the first stage and gradually decreases 
for the successive stages. Such work-done factors can be introduced 
in program AXC03, but for the present compressor they will be taken as 
unity because of the small number of stages and the large blade 
heights. For the same reasons the blockage factor is assumed to be 
equal to unity also. 

The geometry of the blades could be obtained from the information 
of Ref. 5. Values at the mean radius were used for the analysis. 

Bowen (5) determined that the average flow angle at the mean radius 
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after the inlet guide vanes ALE was 20.0 for the design inlet guide 
vane stagger angle of 11.4^, independent of flow rate. From Fig. 12 

it can be recognized that the inlet guide vanes were set at a stagger 

o o 

angle of 10.9 for the tests. Therefore a value of 19.5 was used in 

the analysis program for the average flow angle leaving the inlet 

guide vanes. Reference 5 does not give values of C_ . for the 

blading. It is possible to introduce in the program various values of 

^Dmin rotor and stator blading. Hence it was possible to 

estimate an average apparent value of . for blading by comparing 

Dmin 

measured results to predictions of the analysis program for several 

assumed values of . . 

Dm in 

A prediction of the 3-stage efficiency for different stator 

stagger angles is shown in Fig. 30 for an assumed value of . of 

Dmin 

0.020. The light line connecting the individual peak efficiencies 
represents the calculated operating envelope which could be obtained 
if continuous control of stator stagger angle for maximum efficiency 
were possible. 

For comparison with measured performance a series of calculations 

were made for stator stagger angles of 23.8^, 27.8^, 31.8^, 35.8^, 

39.8^, and 44-3^. For each stagger angle, values of . of 0.000, 

Dmin 

0.006, and 0.008 were assumed. To avoid confusion with the computa- 
tions of actual performance based on measured quantities, which have 
been called runs, each calculation of the analytical prediction 
program has been assigned a two-digit case number. The first digit 
identifies the stator stagger angle of the case as follows: 
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First Digit 


Stator Stagger Angle 


1 


23.8° 


2 


27.8° 


3 


31.8° 


4 


35.8° 


5 


39.8° 


6 


44.3° 


The second digit identifies 
as follows: 


the assumed value of . of the case 

Dm in 


Second Digit 


Assumed Value of C . 

Dmin 


1 


0.000 


2 


0.006 


3 


0.008 



For example, case 42 is the predicted performance at a stator stagger 
angle of 35.8^ with a value of of 0.006. The summary output of 
each case is presented in Appendix D. 



49 



SECTION 6 



DISCUSSION OF RESULTS 

Experimental data for performance measurements were taken 
during 55 hours of running time. Another 81 hours of operation were 
used for miscellaneous calibration and testing, including 44 hours 
for flow rate calibrations. 

The rotor stagger angle was set at 43.8^ for all runs. The 
design stator stagger angle is 28o8^. For the performance tests, the 
stator stagger angle was set at 23.8^, 27.8°, 31.8°, 35.8°, 39.8°, 
and 44.3°. Measured performance is presented in comparison with the 
predicted results of AXC03,and also by establishing dimensionless 
performance parameters. Figures 31 through 36 are graphs of measured 
performance in comparison with the prediction by AXC03 for each stator 
blade angle setting. Figures 37 and 38 are summary plots of maximum 
efficiency and pressure ratio for all stator blade angle settings. 
Figures 39 through 44 are graphs of the non-dimensional parameters at 
each stator stagger angle. Figure 45 is a summary plot of maximum 
efficiencies and associated deflection coefficients for all stator 
blade angle settings. Figure 46 is a summary plot of maximum pressure 
coefficients for all stator blade angle settings. 

On each of Fig. 31 through Fig. 36 the 3-stage efficiency reaches 
a peak value for each stator stagger angle. The referred flow rate at 
which this peak occurs decreases as the stator stagger angle increases. 

Also plotted on each of Fig. 31 through Fig. 34 are the results 
predicted by AXC03 for the three estimated values of . of 0.000, 
0.006 and 0.008. The predicted 3-stage efficiency reaches a peak for 
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each . . 
Drain 



The referred flow rate associated with predicted peak 



efficiency decreases with increase in stator stagger angle. The 
referred flow rate at which the predicted peak efficiency occurs 
decreases with increase in at a given stator stagger angle. 

AXC03 predicts the condition of surge for each case according 
to the criteria of Howell shown in Fig. 29. The surge condition is 
assuraed to exist whenever the pararaeter iq) / € exceeds the 

raaxiraura positive value SMAX which has been inserted in the program. 
The incidence angle for the rotor blade rows becomes larger as the 
through-flow velocity decreases until a surge condition is indicated. 
AXC03 predicts surge at the same referred flow rate for each setting 
of stator stagger angle. This is due to the fixed angle of the 
inlet guide vanes. The assumption has been made that the air flow 
angle leaving a blade row does not change with variations in flow 
rate in the incompressible flow regime. This result was verified by 
Bowen (5). The inlet guide vanes were not rotated as the stator 
stagger angle was varied. At all stator stagger angle settings the 
first stage rotor had the same incidence angle at particular flow 
rates. In the stage-by-stage analysis, "surge in rotor 1" was 
indicated at the same flow rate in each run, regardless of the stator 
stagger angle setting of subsequent stators. The region of referred 
flow rates smaller than that at which surge is indicated by the 
program is marked with a dashed line on Figs. 31 through 34, and is 



less than SMI, the negative value of the parameter which has been 
inserted in AXC03, the program stops the computation and prints 



less than 




When the value of the parameter fp) /^fo determined to be 
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"minimum pressure in blade row xx." The referred flow rate at which 
computation is stopped decreases as stator stagger angle is increased. 
At stator stagger angles of 31.8^ and lower (Figs. 31 through 33) the 
referred flow rate at which computation is stopped is larger than the 
maximum flow rate of the graph. At a stator stagger angle of 35.8^ 
(Fig. 34), the computation stops at a referred flow rate of 45.5 
Ibm psia. At stator stagger angles of 39.8^ and 44.3^ (Figs. 

35 and 36), the program did not establish useful data. At the high 
stator angles the referred flow rate at which computation was stopped 
was lower than the referred flow rate for surge in the first stage 
rotor. The latter flow rate remained unchanged at a fixed inlet guide 
vane angle. 

This difficulty illustrates some of the limitations of program 
AXC03 when applied to an actual machine. The program assumes that 
the compressor is designed properly. Changes from stage to stage 
must be gradual. For example the annulus area may change gradually 
from stage to stage with no adverse effects. The checks for minimum 
pressure and surge limits are applied in each blade row. Surge 
in any blade row is interpreted as surge in the machine, and similar 
limitations are imposed for minimum pressure. The program cannot 
cope with discontinuities similar to those that occur by leaving 
the angle of the inlet guide vanes unchanged. It is possible to 
imagine a situation where the angle setting of the inlet guide vanes 
might have to remain unchanged; for instance, if the inlet guide 
vanes must support the front bearing of the rotor shaft. In such a 
design the first rotor would have higher blade efficiencies than 
succeeding blade rows at large flow rates^ while at lower flow rates 
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the blade efficiencies of the first rotor would be smaller than for 
succeeding blade rows. The 3-stage efficiency of an actual machine 
is an overall parameter which includes the effect of the different 
blading efficiencies for different blade rows. In an actual machine 
the flow will adjust itself to these conditions ^whereas program AXC03 
cannot cope with these peculiar circumstances. For these reasons the 
program did not produce performance data for the higher stator stagger 
angles of 39.8^ and 44.3^. 

The detection of the actual surge point during the tests is 
difficult with the available instrumentation and must be based on 
acoustic phenomena. Incipient surge , probably due to rotating stall, 
was associated with an unmistakable oscillating change in sound level 
produced by the compressor. This effect appeared suddenly even when 
flow rate was decreased slowly. To eliminate the sound due to this 
surge condition it was necessary to increase the flow rate by about 
5 to 10 per cent. A flow rate slightly above the surge point could 
then be reached by gradual throttling. For each run one test point 
was taken in the surge region and one test point as close as possible 
to the surge point. Efficiency at test points in the surge region was 
below 0.90 for all stator stagger angles. The efficiency at the data 
point close to stall rose from 0.89 to 0.94 as the stator stagger 
angle was increased from 23.8^ to 31.8^, and efficiency decreased to 
0.89 with further increase in stator stagger angle to 44.3°. As the 
surge condition is approached the rate of decrease of efficiency is 
set by the relative location on the efficiency versus referred flow 
rate graph of the two data points in the surge region and close to 
surge. At a stator stagger angle of 31.8° which produced the maximum 
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3-stage efficiency and at 35.8°, efficiency decreased most rapidly as 
the surge condition was approached. For other stator stagger angles 
the efficiency decreased less rapidly. It is not possible to make a 
meaningful comparison between measured surge point data and the surge 
condition calculated by AXC03 because surge was predicted to occur at 
the same referred flow rate for each case. 

It is possible to compare measured values of maximum 3-stage 
efficiency with the maxima calculated by program AXC03. For all 
four stagger angles at which computations were made, the maximum 
3-stage efficiencies were predicted to occur at flow rates about 
5 per cent greater than the flow rates at which the measured maxima 
occurred. According to Howell the point of maximum efficiency 
occurs for a value of the parameter )/^/o ^bout equal to 

+0.19 (Fig. 29). At greater positive values of the parameter the 
losses are believed to increase rapidly. Since the measured maximum 
efficiency occurred at a flow rate which was about 5 per cent less 
than was predicted, an analysis of off -design performance of a com- 
pressor stage was made. The method of Vavra (Ref. 1) assumes that 
the air flow angles leaving a blade row do not change with flow rate 
in incompressible flow. The velocity triangle changes caused by a 
reduction from design flow rate are shown in Fig. 47. The subscript 
d refers to the design condition, and the primed quantities refer to 
off-design. The analysis was conducted at the maximum efficiency 
point for a stator stagger angle of 27.8°. The following quantities 
were measured: 

Design flow coefficient = 0.517 

Inlet flow angle cK ^ (fixed) = 20.0° 
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The off-design flow rate was fixed at 5 per cent less than the design 



flow rate, to correspond to the flow rate at which the measured maximum 
efficiencies occurred. The off-design flow coefficient is 
4^^ = (0.95) (0.517) 

The peripheral components of velocity are: 



iL: 

U 



tarx cx I 



= 0.?5* tan o<, 

d 



V^d_ _ U-Vo-Id _ 



u - — u ' - 4^ ”<■ 



The relative inlet flow angles for both conditions are: 

. - ....... Wu,j/U 

r Id 









r. arcian 






- -anc-tan - Cretan 

0.95 <pj 



r j , 

<+, ti.h o(, 

\l 

— _ o<, 

0.9? 4‘j 



With the same rotor stagger angle, the difference in incidence 
angle C between design flow rate and off-design flow rate is the 
difference in 

= 9.'- f.d = ^’'5° = + '-7° 

The difference between the parameter ( L~ L )/^fc design 

flow rate and the parameter at a flow rate 5 per cent less than 
design is 






\o 



+ h7 
I I • i 



= +• 0.15 



where was calculated by AXC03. The measured maximum efficiencies 

occurred at a value of the parameter which was greater 

by 0.15 than the value of the parameter for which the program AXC03 
predicted maximum efficiencies. 

Lieblein*s data were obtained from compilations and reductions of 
cascade test data primarily for blade profiles which were NACA 
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65 (A^q) “ series airfoils built up on equivalent circular arc 

camber lines. The free- vortex blades installed in the compressor 

have a thickness distribution which was derived from the theory of 

Ref. 1 applied to a parabolic mean camber line. If the actual 

incidence for minimum profile loss for the installed blades were 

smaller by 1.7° than calculated by subroutine THEORY, all values of 

the parameter calculated in subroutine CASCAD would be 

larger by 0.15 than they actually are. The flow rate at which the 

maximum efficiency is calculated would be 5 per cent greater than 

the flow rate at which maximum efficiency was actually measured. 

Howell presented results of cascade tests of British C-1 and 

C-2 airfoil shapes for several values of solidity. The C^/C^ . 

D Dmin 

curve of Fig. 29 was obtained by interpolation for a solidity of 
0.8. If the upper limit of the parameter for minimum 

loss were extended by 0.15 from 0.19 to 0.34, the flow rate at which 
maximum efficiency was calculated analytically would coincide with 
the flow rate for measured maximum efficiency. 

The measured efficiency curves are more peaked than the curve 
computed by AXC03. If the curve of Fig. 29 had a greater 

slope in regions where approaches a value of 2.0, the pre- 

dicted efficiency curves would be less flat and would more nearly 
approximate the curves of measured performance. 

It was one object of this thesis to attempt to determine the 

value of C_ . for the installed blade profile. Vavra (1), p. 377, 
Dmin 

indicates that these values usually do not exceed 0.008. For each 

stator stagger angle, analytical calculations were made for estimated 

C_ , of 0.000, 0.006 and 0.008. The results for the smaller stator 
Dmin 
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stagger angles are plotted on Fig. 31 through Fig. 34. The value of 

Ctn j estimated for each stator stagger angle: 

Dm in 

04 - 04 . A 1 Estimated , 

Stator Stagger Angle Dmin 



23.8 


0.008 


27.8 


0.004 


31.8 


0.003 


35.8 


0.007 



The arithmetic average of the estimated values is 0.0055, but there is 
is obviously considerable scatter in the experimental data. 

A summary plot^ of maximum peak efficiencies is presented in Fig. 

37. If continuous control of stator stagger angle were possible, as in 
the proposed SST engine. Fig. 37 would be an operating curve for max- 
imizing efficiency at any flow rate.. The peak maximum 3-stage efficien- 
cy was 0.956, measured at a stator stagger angle of 31.8. Also shown is 

a flatter curve of predicted peak 3-stage efficiencies using . of 

Dm in 

0.006 for the four lower stator stagger angles. The maximum predicted 
value is 0.942 at a stator stagger angle of 30^. The predicted values 
occur at flow rates 5 per cent higher than corresponding measured val- 
ues. The measured values decline more steeply at stagger angles other 
than the angles for maximum efficiency than do the analytically calcu- 
lated valueso The curve relating stator stagger angle to performance 
shows that, in the region near the maximum efficiency, relatively small 
changes in angle are required for large variations in flow rate. 

Between stator stagger angle settings of 27.8^ and 35.8°, a 10 per cent 
variation in flow rate requires only an 8° angle change. In the 
regions more distant from the maximum efficiency point, large angle 
changes are required to handle small flow rate changes. 

Figures 31 through 36 also contain plots of 3-stage pressure ratio 
versus referred flow rate for each stator stagger angle. 
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As expected, the measured maximum pressure ratio occurred at a 
referred flow rate between that for maximum efficiency and for surge. 
It was not possible to determine maxima for the calculated data of 
AXC03 because for each case the largest pressure ratio occurred at 
the flow rate tagged "surge.** The program does not compute perform- 
ance at lower flow rates. The check for surge occurs in subroutine 
CASCAD. Based on cascade data^Howell stated that unacceptable losses 
occur at values of greater than 1.25. This conservative 

criterion is particularly applicable in a design where a wide range 
of flow rates is not required. However for the purpose of evaluating 
pressure ratio for an existing machine, the criterion appears to be 
too restrictive. It would be more advantageous if AXC03 used the 
criterion to predict surge but then proceeded to complete the calcu- 
lations at that flow rate rather than stopping computation* 

The measured pressure ratio was in general by 0.25 per cent 
greater than predicted by AXC03, and the measured ^ about 

11 per cent greater. Total pressure is relatively simple to measure. 
The yaw and pitch errors of the 3-hole probe are very small. Any 
error in total pressure measurement is expected to produce a total 
pressure lower than expected. Since the predictions are lower at 
every flow rate, it is possible that the difference in blading 
between the installed blade profiles and the NACA airfoils, which 
are represented by Liebleins's data, may be the cause of the differ- 
ence as was suggested in the discussion of 3-stage efficiency. A 
summary plot of maximum 3-stage pressure ratio on an expanded scale 
is presented in Fig. 38. The peak value of maximum pressure ratio 
occurred at a stator stagger angle of 31.8^ 
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While writing the thesis, an error was discovered in Eq. 1, which 

converts the reading of the mercury barometer P in in. of mercury 

Bar 

to in psfa and corrects for the temperature variation of the 
A 

specific gravity of the liquid. The incorrect equation 

= Pg^^ ( 0 . 492.5 -t- 0.5/10000)144 (psfa) 

was used. At 70.0°F the constant by which was multiplied was 

71.457 psfa/in. of Hg. Reference 10 lists the specific gravity of 

00 

mercury sy for the temperature range 0 F to 150 F as 

sy = /5,(>38 “ '^5'4 X ) (g/cm^) 

The conversion factor is 

Ct - (psfa/in. of Hg) 

At 70^F the correct conversion factor is 69.55, or the atmospheric 
pressure on which all calculations were based was 3 per cent higher 
than the actual atmospheric pressure. By Eqs. 3, 4 and 7 the cal- 
culated static pressure P^ and air density at the permanent 
Pitot-static tube, and the total pressure after the inlet guide vanes 
Pt^gafWere 3 per cent higher than actual also. By Eqs. 6, 8 and 10 
the calculated velocity V^, volume flow rate y and referred 

flow rate w^ were 0.983 of the actual values. The calculated pres- 
sure ratio P « was 0.998 of the actual value by Eq. 11. By Eqs. 9 
r3-s 

and 14 the values of weight flow rate w and 3-stage efficiency 
were calculated to be 1.2 per cent higher than the actual values. 

The test results are presented in non-dimensional form in 
Figs. 39 through 44. Plotted against the flow coefficient , the 
efficiencies follow the same pattern as when plotted against referred 
flow rate, since EqSc 10 and 13 differ only by a constant. Figure 45 
is a summary plot of the maximum 3-stage efficiciency for each stator 
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stagger angle versus flow coefficient and is similar to Fig. 37 for 
the same reason. Figure 46 is a summary plot of maximum pressure rise 
coefficient for each stator stagger angle. The pressure rise reaches 
a maximum at a blade angle of 33*^ and a flow coefficient of 0.475. The 
curve of corresponding stator stagger angles has an almost constant slope. 

Theoretically^ the deflection coefficient ^ 0 ^ represents the 
change in peripheral flow components in a blade row divided by the 
peripheral speed. Figure 47 is a stage velocity triangle showing 
design and off-design conditions in a compressor stage. At the design 
condition 



r 









d 



u u 

An assumption in this analysis is that the air leaving angle from a 
blade row does not vary with incidence. Hence the angles j and 
remain constant. The relation between the off-design deflection 
coefficient L and the off-design flow coefficient ^ is fixed by 



similar triangles 
^ * 

I — L 



<p 



I - i 






and 



T' = / - ( I - TJ 









or L is expected to vary linearly with . 

At flow coefficients greater than that for maximum 3-stage 
efficiency^ calculated deflection coefficients varied linearly 

with flow coefficient in Figs. 39 through 44. For most runs the 
graph of Z ay broke sharply from the linear segment at the maximum 
efficiency point. For each case the curve continued at a reduced 



60 



slope in the region between the break and prior to surge. Experience 
seems to indicate this identifies a region of incipient stall. As 
the blades become fully stalled, the nature of the flow changes 



drastically. It has been suggested that this change is reflected in 
a sharp increase in the slope of the deflection curve. The author 
was reluctant to operate the compressor in the region of stall for 
extended periods of time because of related vibration and noise prob- 
lems. However a slight upturn of the deflection coefficient curve 
was apparent at stator stagger angles of 23.8^ and 27.8^ (Figs. 39 



will be necessary to make an approximation for the theoretical power 



and 40) . 



The deflection coefficient X is a measure of work input. It 



, . Equation 11 can be rewritten 
Wth 






(psfa) 



Making a first order expansion, Eq. 13 can be rewritten 




(ft-lbf)/sec 




(f t-lbf ) /sec 



Substituting P,, , in Eq. 12, the efficiency becomes 
wth ^ 





Rewriting Eq. 15 




(f t^/sec) 
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and introducing Eq. 16 into Eq. 17, there is 



or the deflection coefficient is a direct function of the torque 

T or work input* 
rq 

Figure 45 is a graph of the which was calculated at max- 

imum 3-stage efficiency for each stator stagger angle. The minimum 
of these deflection coefficients occurs at the maximum efficiency. 
At maximum efficiency the minimum work input was required. 
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SECTION 7 



CONCLUSIONS AND RECOMMENDATIONS 



It has been shown that raaxirnuin efficiency can be achieved for 
each flow rate. The computer program AXC03 can be used to predict 
efficiency within 2 per cent^with an error of 5 per cent in flow 
rate. Pressure ratio may be predicted within 0.25 per cent. 

It is recommended that further measurements be made at increased 
rotor stagger angles. Inlet guide vane angles should be varied to 
simulate the effects of preceding blade rows. Blade shapes with 
higher blade loadings should be tested. 

Scatter in the data may be caused by variations in inlet flow 
and accentuated by oscillations of the pressure gauge. A proposed 
improved inlet duct is described in Appendix E. 
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FIGURE 1 



COMPRESSOR INSTALLATION 
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FIGURE 9 COMPRESSOR CASING 




FIGURE 4 ROTOR BLADES, DRUM AND SHAFT MOUNTED IN 

SUPPORTING STRUT ASSMBLIES 







FIGURF 5 



ASSEMBLY OF ROTOR IN LOWER CASING 
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FIGURE 6 COMPRESSOR SECTI6NAL VIEW 
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FIGURE 7 LOCATION OF RADIAL SURVEY HOLES 
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FIGURE 8 DETAIL OF FREE -VORTEX ROTOR 
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FIGURE » DETAIL OF FREE-VORTEX STATOR 
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free-vortex stator blade 



c-/ r/xa> sscro^ 
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FIGURE II ROTOR STAGGER ANGLE ADJUSTMENT 




FIGURE 12 EXTERNAL STAGGER ANGLE ADJUSTMENT 
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FIGURE 13 INLET BELLMOUTH 



7b 




i 



FIGURE 14 INLET PITOT- STATIC TRAVERSE 
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FIGURE 15 DETAIL OF INLET PITOT-STATIC TRAVERSE 





FIGURE 16 THREE-HOLE PROBE DETAIL 
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GURE18 PERMANENT PITOT -TUBE AND THhtt 
WOLE PROBE AHEAD OF FIRST ROTOR 





FIGURE 19 PRESSURE READOUT INSTRUMENTATION 
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FIGURE 22 TRAVERSE CARRIAGE AND THREE — 
HOLE PROBE BEHIND THIRD STATOR 
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FIGURE 23 RELATION OF TRA/ERSE CARRIAAE TO SURVEY HOLES 




FIGURE 24 TORQUE METER DURING CALIBRATION 
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FIGURE 25 TRAVERSE OF INLET' DUCT RUN D-l 
VOLUME FLOW RATE 512 FT^/SEC 
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FIGURE 26 TRAVERSE OF INLET DUCT RUN 0-7 
VOLUME FLOW RATE 310 FT^/SEC. 
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FIGURE 28 TEMPERATURE-ENTROPY DIAGRAM OF 
COMPRESSION PROCESS IN 3-STAGE. 
MACHINE 
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FIGURE 3( 3-STAOE PRESSURE 
RATIO a 3-STAGE EFFICIENCY 
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APPENDIX A 



CALIBRATION OF INSTRUMENTS 



The measuring devices were calibrated to determine their accuracy 
and precision. 

The Berkeley digital counter is located near a Hewlett-Packard 
digital counter in the control room» At a given flow rate, there 
was no difference between the two instruments. The Hewlett-Packard 
device was recently calibrated with a cesium wave length frequency 
standard in the Standards Laboratory of the Electrical Engineering 
Department of the Naval Postgraduate School. The speed measured is 
considered accurate to 0.5 rpm. 

The Baldwin-Lima-Hamilton torque meter was calibrated statically 
as illustrated in Fig. 24. A bar was attached to the drive end of 
the shaft and was clamped to the stand. A symmetrical lever bar was 
attached to the rotor end of the shaft. The lever arm is 20.0 in. 

Two weight pans were adjusted to 1.44 lb. each with lead shot in a 
plastic bag. The weights themselves were checked on the Toledo 
Precision Scales of the Cascade Laboratory. On the first run, dif- 
ferent readings were obtained while loading and unloading the pans. 
This problem was overcome in subsequent runs by tapping the stand with 
a mallet prior to recording the reading. This simulated the vibration 
of dynamic operation which eliminated any frictional effects. Good 
agreement was obtained during tests where the loading was increased 
and decreased. The possibility of temperature effects on the strain 
gauges was explored in a third test. The torque meter cannot be 
calibrated while the motor is running, but it was calibrated 



113 



immediately after shutdown. The shaft was only slightly warmer to 
the touch than ambient temperature. 

The data of the three calibration runs are shown in Table A-1. 
The results are plotted in Fig. A-1. Run A- 3 revealed no temperature 
effect. A fourth run was made later for checking purposes and did 
not reveal any discrepancies. The calibration tests established a 
constant for the torque meter of 4.1565 foot-pounds of torque per 
pound reading of the meter. The maximum relative error among the 
runs is 0.6 per cent. The meter may be read to within 0.03 pounds, 
or 1 per cent of the meter reading. 

The newly acquired low range bourdon tube was calibrated also. 
Run B-1 was made by using the mercury manometer board in the control 
room of the Compressor Laboratory. The transsonic turbine test rig 
was used to provide the pressure difference for the higher pressure 
range. The present compressor tests required pressure measurements 
down to 0.6 in. of water. Runs B-3 and B-5 were made at low pres- 
sure difference with a Merriam water micro-manometer. Pressure was 
supplied from a static source. 

The data of the three calibration runs are shown in Table A-2. 

The results are plotted in Fig. A-2. The calculation of velocities 

2 

was set up using pressure difference in Ib/ft . Measured pressures 

are therefore very small numbers. An inverse calibration constant 

was used. The bourdon tube has a constant of 240.423 counts per 
2 

lb/ ft . The maximum relative error among the calibration runs was 
0.12 per cent. The meter may be read to within 5 counts. The 
meter’s servo system tends to overcorrect. In the manual null mode 
the meter needle tends to wander a similar amount. The five count 
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reading error amounts to 1 per cent of the velocity head reading in 
the inlet duct during flow rate calibration. The error is 0.5 per 
cent of the velocity head reading at the permanent Pitot tube. The 
error is only 0.05 per cent of the measured pressure difference 
across the three stages. 

The Pitot-static tube used to traverse the inlet duct during 
flow rate calibration is a modified Prandtl-type tube. Another 
Pitot-static tube was attached to the traverse shaft and the tubes 
compared. Sufficient difference existed between the two to justify 
further tests. The low speed calibration tunnel in the Cascade 
Laboratory was modified to receive the Pitot-static tube from the 
inlet duct. 

The data of the calibration runs are ffhown in Table A-3. The 
results are plotted in Fig. A-3. Each point represents the average 
of four readings. Run C-3 indicated that the velocity head of the 
Pitot-static tube from the inlet duct was at most 1.1 per cent less 
than that of the calibration tunnel. Run C~4 indicated that the 
static pressure below atmospheric was at most 1.8 per cent less than 
that of the calibration tunnel. These maximum errors occurred at the 
lower pressure differences. It is probable that the effect recorded 
is in part the reading error due to bourdon tube oscillation. 
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TABLE A-1 




t i 

!• 



1 



CALIBRATION OF BLH TORQUE METER 



RUN A 


-2 




RUN 


A-3 


APPLIED TORQUE ! 
FT-LB i 


METER READING 
LB 


APPLIED TORQUE 
1 FT-LB 


j METER READING 
1 LB 


i 

4.167 


0.995 


1 19.792 




4.745 ! 

1 ■ V 


8.333 ' 


2.005 


’ 39.583 




9.620 ij 

'! 


12.500 


3.000 


59.375 




14.540 i! 

i 1 


16.667 


3.985 


77.708 




18.445 '! 

\ 


20.833 


4.995 


92.708 




22.400 H 

’ . i 


' 25.000 


5.985 


( 

109.375 




26.345 1 

1 


1 

29.167 

i 


6.980 i 


i 126.042 




* 1 • 
30.285 1 


' 33.333 


7.960 * 






1 


i 37.500 


8.940 


1 

i 




I, 


41.667 


9.970 


! 

1 




J f 


i ” 






45.833 


10.925 


( 


RUN 


A-4 1 


50.000 


11.950 


i APPLIED TORQUE 
FT-LB 


! METCR READING ' \ 

1 , 


54.167 


12.915 ; 


19.792 




4.77 j 

< 


58.333 


13.930 


; 38.958 




9.610 i 

1 


62.500 


14.905 


59.792 




14.545 


^ 66.667 


15.925 


i 76.458 




18.510 


70.833 


16.935 i 


j 93.125 




• 22.435 

j i 


75.000 


17.945 ! 

1 


1 

! 109.792 




j 26.405 

1 


79.167 

1 


18.910 1 

! 


( 

126.458 




! 30.370 ' 

i 1 


83.333 


1 

19.970 

, 


130.675 




' 31.340 

j 

1 
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Meter Reading^ Pounds 



FIGURE A-l TORQUE METER CALIBRATION 






Average Colibrotion Conetont 4.I56S FT-LB Per Pound 
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TABLE A- 2 



CALIBRATION OF 0-12 in. Hg BOURDON TUBE 



RUN B-1 



PRESSURE 
INCHES OF 



METER 

READING 



RUN B-3 



PRESSURE 



METER 



RUN 



PRESSURE 



INCHES OF i READING ' INCHES OF 



MERCURY 


COUNTS 


' ’ WATER i 


COUNTS ; ; 

; j 


WATER 


10.260 


174,300 


' 13.978 


17.409 i 


14.519 


8.290 


139,450 


12.459 


15,518 


12.914 


7.445 ! 

1 


125,750; 


11.062 


13,784 ,j 

1 1 


11.397 


6.065 1 


102,900 


9.853 


12,276 1 


9.913 


5.235 


88,200 


6.993 


8,737 ; 


8.471 


4 . 440 


74,800 


1 6.198 


7,738 


6.989 


4.035 


67,950 i 

1 


4.992 

1 1 


6,230 1 


5.471 


3.320 


! 57,100' 


! 3.996 

t 


4,993 1 


4.493 


2.415 ' 


' 40,900 


i 

i 3.036 

j 


3.794 1 


3.470 


1.000 


16,300 

1 


! 2.195 

j 


1 

2,741 1 


2.401 




1 

1 


!l 1.602 


1,997 


1.512 






0.785 


982 


0.620 



B-5 

I METER 
READING 
i COUNTS 

i- 

! 18,218 
i 16,119 
14,222 
12,368 
10,560 
8,715 
6.826 
5,613 

I 4,398 

i 

I 2,994 
1,885 
762 
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FIGURE A-2 



BOURDON TUBE CALIBRATION 






Avtrage Calibration Constant 240.423 
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TABLE A- 3 



\ 



CALIBRATION OF INLET PITOT-STATIC TUBE 





RUN C-3 




RUN C-4 


TOTAL PRESSURE 




STATIC PRESSURE 


MINUS 


STATIC PRESSURE 




BELOW . 


ATMOSPHERIC 


PROBE 


i 

i TUNNEL 


PROBE 




TUNNEL 


2486.25 


2508.75 


2450.0 




2455.0 


2157.5 


2182.5 


1 2190.0 

! 




2190.0 


1413.75 


1421.25 


) 

! 1422.5 

j 




1425.0 1 


1002.5 


1011.25 


1 

1 1087.5 

5 


) 

1 

1 

1 


; 1012.5 ! 

1 1 


793.75 


797.5 


1 

1 795.0 


) 

1 


810.0 ' 




i 

» 


f 

» 

f 

j 

i 


i 

1 

j 

i 

i 


1 

i 
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Rrobe Static Pressure Below 

Atmosphere, Counts Probe Velocity Head , Counts 



FIGURE A-3 INLET TRAVERSE PITOT-STATIC CALIBRATION 





Below Atmosphere, Counts 
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APPENDIX B 



FLOW RATE CALIBRATION DETAILS 



Particulars of the flow rate calibration scheme which may be of 
special interest to a restricted number of readers are included in 
this Appendix. Section B-1 gives details of the construction and 
testing of the variable length capillary damping devices. Section B-2 
gives a listing of the program ONRFLO, and Section B-3 presents sample 
output of program ONRFLO. 

B-1 Variable Length Capillary Damping Devices 

The problem of oscillating pressures at the traversing Pitot- 
static tube in the inlet has been described. Pressure oscillations 
in the tubing were damped with alternating sections of smaller and 
larger diameter flow passages. Figure B-1 is a picture of such a 
device. The device consists of three lengths of needle stainless 
steel tubing, of 0.020 in. inside diameter, which are 4 ino , 6 in., 
and 8 in. long. They are connected as shown to polyethylene tubing 
by brass **T*s.** By clamping the tubing at various positions it is 
possible to obtain total capillary lengths of 4 in. , 6 in. , 8 in. , or 
18 in. Figure B-2 is a schematic showing the operation of the device. 

Two variable length damping devices were used, one each in the 
total and static pressure lines from the traversing Pitot-static tube 
in the inlet. Two effects of these devices were observed on velocity 
head, the difference between total and static pressure. As the length 
of capillary tube was increased, the fluctuation of the pressure gauge 
was decreased. But as the length of capillary tubing was increased. 
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the time to come to a steady reading was also increased. Table B-1 
presents the results of tests of various lengths of capillary tubing. 

Immediately before taking the velocity head reading the traverse 
probe was placed at a new radial position. This task took about 20 
sec. to perform. The best precision for a time to damp of 20 sec 
or less was used. The combination of a total pressure capillary 
length of 18 in. and a static pressure length of 8 in. satisfied this 
requirement. 

Later in the investigation, when the simplified pressure selec- 
tion system was used, it became possible to use both 18 in. lengths 
of tubing and have a time to damp to a steady reading of less than 
20 sec. It appears that the 12-channel pressure scanners increased 
the time required for the system to stabilize. This was probably due 
to admittance of air from other channels during sequencing of the 
scanner. However as the survey progressed, flow rate measurements 
were made at lower flow rates. As the machine was operated nearer to 
the stalled condition, greater pressure oscillations were observed. 
The meter precision was no better than + 15 counts, and sometimes as 
much as + 50 counts. 

Program QNRFLO 

The flow rate calculations described in Section 3 were accom- 
plished by program ONRFLO which is included as Table B-2. Since the 
calculations are straightforward and have been described earlier, 
they will not be repeated here. 
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The input format is as follows: 



Card 1: 


NCASE 


run number (arbitrary) 


(I7,3F7.2) 


PEAR 


reading of barometer (in.) 




TEAR 


temperature at barometer (°F) 




TO 


temperature at inlet screen (°F) 


Cards 2, 3, and 


4: R 


the 23 radii at which data were 


(8F1-0.2) 




taken (17o75, 17. 5 0, . . .4.00, 0.00, 
4.00. .ol7. 50, 17. 75) (in.) 


Cards 5 through 


50: QI 


velocity head in inlet duct (counts) 


(4F7.0) 


PSI 


static pressure with reference to 
atmospheric in inlet duct (counts) 




QO 


velocity head at permanent Pitot- 
static tube (counts) 




PSO 


static pressure with reference to 
atmospheric at permanent Pitot- 
static tube (counts) 


The radii of cards 2, 3 and 4 


are repeated internally to account for 


the cross traverse. Cards 5 


through 27 introduce measurements made 



at the data points for one traverse, and data for the cross traverse 
are introduced on cards 28 through 50o 

B-3 Output of ONRFLQ 

Table B-3 is a sample of the output of ONRFLO. It is for the 
final calibration, run 14. The first data block is a list of input 
data. For each point on the two traverses, QI , PSI, QO, and PSO are 
listed. The second data block is a list of computed velocities VI, 
VO, and VIAD. The third block shows how the **no boundary layer** 
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flow rate, VFT, is computed by summing the product of VA and A. 

The fourth block shows how the flow rate deficiency, VFD, is computed. 
It lists the inputs to and outputs from SUBROUTINE QTFE, which per- 
forms the integration by the trapezoidal rule. The final statements 
list the final computed flow rate VFC and the final calibration 
constant CCC. 
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FIGURE B-l VARIABLE LENGTH CAPILLARY DAMPING DEVICE 



I 2 



VARIABLE LENGTH CAPILLARY 
DAMPING DEVICE 
OPERATION 




MATERIALS 

CAPILLARY: NEEDLE STAINLESS TUBE 0.020" LD. 

valve: polyethylene TUBE 8 LAB CLAMPS. 



CAP. 

LENGTH 

DESIRED 


VALVES 


A 


B 


C 


D 


4" 


CLOSED 


CLOSED 


OPEN 


OPEN 


6" 


OPEN 


CLOSED 


CLOSED 


OPEN 


8" 


OPEN 


OPEN 


CLOSED 


CLOSED 


18" 


CLOSED 


OPEN 


OPEN 


CLOSED 



FIG. B-2 



127 



PRECISION a TIME TO DAMP FOR 
VARIOUS LENGTHS OF CAPILLARY 
TUBE IN TOTAL AND STATIC PRESSURE 
LINES. 



CAR LENGTH IN. 


PRECISION, 


TIME TO 
DAMP TO 


TOTAL 


STATIC 


+ COUNTS 


STEADY. 

SEC. 


18 


18 


1 


45 


18 


8 


5 


20 


18 


6 


8 


10 


18 


4 


10 


5 


8 


18 


5 


45 


8 


8 


5 


35 


8 


6 


10 


20 


8 


4 


10 


5 


6 


18 


5 


30 


6 


8 


10 


15 


6 


6 


10 


10 


6 


4 


20 


10 


4 


18 


15 


30 


4 


8 


25 


15 


4 


6 


30 


10 


4 


4 


40 


8 



TABLE Bl 
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Table B-2. Listing of Program OKRFLO (cont.) 



00 



< 



< 



X 

< 



00 

h- 

U2 

H-O 

^/^oOLf^lr^u^lr^u^ou^moou^ou^oooi^uf^u^ool^kOOOOu^ou^^/^u^L^^olJf^Ol/^Ol/M/M/^lrM/^u^ 
LU ^ LH >t ^ ^ '4' ^ U^ sf '4’ >t 't' ^ ^ '■t '4* ^ ^ 

UJ QC (NJ r^J (N4 (NJ r^J (NJ <NJ CSJ Csj (NJ e\J (NJ (N fNJ (\J <N) ru (N (NJ <\J <N CNJ (NJ <\J <\J (NJ fsj (M (M rsj <\J (\J <\i rvi fvi rsj (\J rg ro rvj fSj fSJ <NJ (Ni (NJ fNJ 
CDZ) 

QOO 
Qi 00 
CLLU 
QC 

a 



2 • H-00 

>0 5: 

• oZ 

Q K CL 03 •••••• 

h- < # -jOuv/M/MrsoooirkLr>oi/>ooootfMr>OLnLr>oiAoir^oir>Lnou^ooooi/>04riirki/>oirycoLryooiry 

< a: oc uJoooc^^•Lf^»oc^^-oocoooo^^voa^^•ooO"^4•oo^^^oooc^a^oocx5^-oo^roo^^oo>4*'Of^QO^-a'ooo^oc^ 

or XI > h-r^>>f^xr^r^^-r^h-N-xr^r^r^r^h-r^oo>h->r^r^h-r^r^h-r^QOoor^r^c30'0'Cr^>r^r^Kcx)oooor^ 

tmmt ^m4 ^m4 pmi4 m 

cn u o ujQ 

< < (\i x< 

^ u ^ Olu 

OCX 

z Lu II a 



« 



3 

a 

a 

gc 

CL 

O 

X 

oc 

3 



c 

3 

X 

o 

00 



:s 

o 



QC 

3 

oO 

00 

UJ 

CL 



O 



00 

^3 

h-O 

<o 

h- ir>ixvoOiryixvirsOirsirvirsLnLnoooipiryoooirii/>ir\iryOir>oooixsOoi^Lnooir>oirvLrMryLrMnLnLr 
00 •‘C\jvOvcr^>t>4mrnh-a^^oro>xvoi^r<^mN4m>t^Xfrio^'rirvvOCN>oofor^fNifornxmm>t 

LU X 00 X X X 00 (X) X X X CT' CO 00 X 00 00 00 X 00 ® r- oc 00 X X 00 X X 00 X X QC CT' 00 X X X X X X 00 X 

h-ar 
LU3 
^00 
zoo 

*-^LU 

QC 

a 



Qt 

3 

3 

CL 



< 

> 

z 



.OO 



ty^ 

I 

ro 



3 



I 



3 

X 

QC 

3 

IX 



OZ 

03 

3 3 iTM^ irv 1 ^ LTy o o m o ufMn o o o o iPk o in o o iPk LrunirMrMrv irs iTMTk o o iroTi o o o o o iTi o o Lf> o in m LA 
Luo>t0voxooocr '\jnj^irk4-^iry<Noa'0>roxc\j^Lnxooo^r\jmfnxrnrom-^xxot7'r^‘4'<M>t 



►-0 

LU< 

30J 

ZX 



5 



< 

3 



3 

CL 

Z 



OO 

^ 1 ) 

x5 j 

OZ * 

QC'- I 

LL LAOir'OLAOinoooo ooooinOLAOLrkOi^LnOLnoinoLrioooo oooomoi^oi^oirk 
LA rg O > ITN c\] OO O O O O O OO fVJ LA O f\j UA A- > m rvj on- iO O O O O OO O O O (N m h- o CNJiXn- 
oOUJ •••••#•••••••••••••••*•••••••••••••••«•••••••• 

3Z>>n-n-v0>C vO vOfMX >r 0 4*XfM vO sOnO xr^n-n^n-r^n-n-n-NO vDsO 0'Njx>r o>4X(MX 

*— •* i^H 0—4 t—4 »..< r—H »— 4 f—^ •— < 

03 
<X 
ocuj 



z 

UJ 

o 



132 



Table B-3. Sample Output of Program OWRFLO- Input for 



CALCULaTFO VtLUUTieS, FEET PER SECOND 



LU 
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00 CM (M (\j m (\j m o fo o >t m CO CO UMn O' O' 
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Table B-4. Sample Output of Program ONRFLO- Velocities for Run lU 



o 

z 

o 

o 

OJ 

00 



Hi 

CL 























o 












o 












u. 






UJ 






''s 


















< 






c 






cat: 






z 








NtO^Ornrn 




o 




UJ 


:x 


r^fsjoooo 




u 




-J 


o 


>tsOfNJOsjO 




LU 




z 


«j 


• » • • • 




00 




1— f 


u. 














ltn 




CC 




z 








UJ 






< 






a 
















OO 


O 












►- 




ru 


LU 




00 








UJ 




0-4 






00 


UL 




X 






• 






UJ 






O' 


O 


0^ 








r- 


»— < 


LU 


O' 


UJ 






CO 


V 


UJ 


h- 






3 


< 


X 


< 






O 




< 


q: 


r^iTNUorvjo 


M 










'-•'OO' 00 






> 






• • • • • 


>- 


O' 


a: 


> 


Q 


^ sO'-^OO' 


1— 


LO 


< 


O' 


-J 


POf^oro 


»— 1 


O' 


c 


< 


UL 


>**^ 


U 




z 


c 






CJ 


• 


3 


z 










3 


3 


o 




UJ 




CQ 


o 


h- 




> 






CQ 


u 








Z 




lU 




UJ 


OO 




a 


00 




o 


1— « 




z 






c 




V 








CC 




u 


(0 






a 


z 


z 


z 








< 


UJ 


•— t 


< 


f^00O>t fO 


OJ 


H- 


t— H 


5: 


UJ 


00 O' O' O' O' 


o 


OO 


o 


3 


a: 


oomor^ 


<r 


z 


1— « 


00 


< 


• f • • • 


QL 


3 


LL 


00 




O O (M 


LU 


O 


LU 


< 






> 




O 








< 


z 




UJ 








3 


15 












o 


< 








1— 




CJi 








< 


LL 










ot: 










'O 


3 


LU 


a 


> 




O 


1^ 


z 




h- 




o 




3 


u. 






• 


< 


-J 




o 




O' 


O 


O 


UJ 


o 




LTN 




> 


z 


u 


rnoocNjcM 


ro 


? 






UJ 


nOto— * in o 




Of 


LL 




> 


osjr^moo 




LU 


3 


o 








> 




> 


UJ 


»— 1 ^ ^3 




< 


Z 




o 


iTMiOLnin 


II 


-J 


a 


U. 


< 










a 






LU 


> 






LL 






q: 


< 


z 


> 




< 


< 




:3 


< 




ut: 


o 


3 


— < 








z 


U 










3 


3 


< 






a 


O 


< 










CO 


O 




o: 




LL 






'O 


a 






o 






H- 




UJ 


z 




< 


o 


^(Nj''nvr o 


Z 






CJ 


LU 




3 








OO 






UJ 










a 


X 










> 


h* 











o 


-4- 

H 








z 










o 


a 








o 

3 


K 








00 


u 










o 








OC 


Ch 


> 






3 




UJO 

h-Z 




Nt 


a 


1 


<UJ 






h- 


+> 


a:i— » 






o 


w 


o 


OOOO-^>t0C00 


• 


o 




3C*-« 


• ••••••• 




3 


vj 
r ^ 


OLL 


OOOOOOO'-^ 


m 




V—' 


3LU 




fO 




G 


U-O 








O 








Q 

Z 


•H 








o 


U 






II 


o 


rO 








3 


•H 






LU 

h— 


00 


o 






< 


Of 


1 






of 


3 


Q 








Q. 


3 






:x 






o 




o 


h- 


§ 


LLZ 




3 


3 


O 


o< 


h- m O' 'O m f<^ 


U. 


3 




CC 


•:'Ooj>tooooLnoo 




3 




OJO 


ororOsorvjr-r-3 


LU 






3LU 




z 


O 




3h— 


033000U>0'0 


3 


1—4 


o 


<Z 


^r\iO 


3 


CQ 


u 






O 


3 


Ph 






> 


O 


Ch 










o 






3 










< 


r^ 


i£ 






Z 


m 








1— * 

UL 


o 


f 








vf 


<D 










3 


> 








& 


xo 




O' 


00 




u»Z 


ro r— 4 m o 


in 


»— 4 


CO 


^LU 


LnLnr>-rn>^>tr*- 


ro 






O— • 


0^,-4 0(jvf^CNJ >t 


• 


I— 




OO 






z 


lf\ 

p4 




00000'~»r0>0 




< 


LLLL 






3 


>LU 






00 


dJ 

3 


O 




II 


z 

3 

O 






> 










O 


z 








z 


3 








3 






»• 




l<-4 


h- 




00 




o 


z 


f 


LU 




*— » 


< 




XX 




3 


CO 




OO 




3 






Of z 




3 


3 




LL— • 


omou^OLrioiri 




< 




»> 


Of\iir>r*-O(\j»jor^ 


Ji: 


O 




00 UJ 


• •••••»• 


a 






3Z 


\0>o f'— r— r"* 


3 


3 




1— 1 


i«-J »«-4 »*-J fo-4 ^ r— 4 


3 


< 




03 






Z 




<aoc 




3 


— « 




TtlLU 




Z 


3 




h- 




3 






z 


134 


3 


3 




3 


O 


X 




o 




> 


3 





APPENDIX C 



PERFORMANCE MEASUREMENTS 

It was decided to take flow measurements at the mean streamline, 
using the radial survey holes. The assumption was made that measure- 
ments made at this point accurately represented the flow in the 
entire flow annulus. The truth of the assumption is verified in 
Section C-1. Section C-2 presents details concerning the data re- 
duction program ONRETA. 

C-1 Verification of Assumptions 

The assumption that measurements made at the mean streamline 
using the radial survey holes accurately represent the flow in the 
entire flow annulus was verified in three ways. The total pressure 
measured at the mean radius was shown to be nearly equal to the 
integrated average of total pressure at each radius. It was demon- 
strated that the total pressure measured at the peripheral location 
of a radial survey hole was nearly equal to the integrated average of 
total pressure at each peripheral location in a blade channel. The 
flow was shown to be nearly axisymmetric. 

Concerning the radial variation of total pressure, the existence 
of a boundary layer at the outer radius in the inlet has already been 
shown. It was further shown that a boundary layer existed ahead of 
the first stage and behind the last stage. With the stator stagger 
angle at 27.8^ and the flow rate at about 400 ft^/sec. , a radial 
traverse was made at the radial survey hole labeled 37^27*30** in 
survey plane 8 of Fig. 7. Data points were established each 0.3 in. 
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from a radius of 11.1 in. to 18.0 in. After recording the barometric 
pressure ambient temperature the following measurements 

were taken: 

- total pressure with reference to atmospheric after third 

stator (counts) 

- static pressure with reference to atmospheric after third 

stator (counts) 

- velocity head after the third stator (counts) 

The total pressure at the inner radius could not be measured but was 
assumed to be the static pressure at 11. 1 in* With the relations 
listed on p. 35 these measuring data establish the following 
quantities: 



p. = ' 


P|i»r 


( 7I.4C.7) 


(psfa) 


(1) 






45^7 


(°R) 


(2) 


Ptj!. ° 


Pa - 


Ptd 

ZA0.4Z3 


(psfa) 


(3) 


Psda' 


Pa - 


Psd 

^4o.4^3 


(psfa) 


(3) 



With the assumption that the total temperature after the third stator 
is nearly equal to the inlet total temperature 




0 . 002.376 5 | 8.7 

CK.G'I 0(144) 



^^^5,S2xlO j(lb sec^/ft^) 

Tto 



(4) 



I 



%d 

240. 4Z3 



(psf) 



(5) 
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^ IJI 



(ft/sec) 



( 6 ) 
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The average total pressure was calculated from the integral expression 

/ 8.0 






P 



fO.Q 






18.0 



IV t 



/o.B 



Both integrations were performed by the trapezoidal rule. Figure C*-l 

is a graphical representation of the integration. The integrated 

average total pressure 2184.20 Ib/ft^ while the total pressure 

2 

measured at the mean streamline was 2184.34 Ib/ft . The relative 

error between the two pressures is 0.006 per cent. A second radial 

survey was made at the radial survey hole labeled 149^57*30" in the 

plane called S. P. 8 on Fig. 7. The integrated average pressure for 

2 

this survey was 2183.92 Ib/ft . The total pressure measured at the 

2 

mean radial position was 2184.32 Ib/ft . The relative error between 
the two pressures was 0.018 per cent. Since the pressure measured 
at the mean radius was nearly equal to the integrated average of pres- 
sure at each radial location, mean streamline values were used to mea- 
sure compressor performance. 

The traverse carriage was used for peripheral surveys. Since it 
would be necessary to remove the first stage rotor to survey in the axial 
plane ahead of the first stage, the investigation was confined to the 
plane behind the third stage stator where measurements could be made simply. 
The locations of the radial survey holes with reference to the indicated 
meridional angle of the carriage are shown in Fig. 23. 
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There are 32 stator blades. Each blade channel is 11.25^ in peripheral 
direction. With the stator stagger angle at 27.8^ a peripheral survey 
was made at the mean radius, 14.4 in. The average reading over the 
entire 15° of carriage travel was 9256 counts. The reading taken at 
the radial survey hole marked 37°27*30" on Fig. 7 was 9246 counts, a 
relative error of 0.11 per cent. The average over the first 11.25° of 
carriage travel, representing one blade channel, was 9218 counts, a rel- 
ative error of 0.30 per cent. Figure C-2 illustrates the results of 
this mean radius survey. A second survey was conducted near the hub 
at a radius of 11.2 in. The average pressure reading of the traverse 
was 9532.47 counts. The pressure measurement at the radial survey hole 
was 9569.23 counts. A relative error of 0.36 per cent existed between 
the two pressures. On the basis of the surveys made behind the third 
stage stator the measurements taken at the radial survey taps were 
accepted as accurately representing the flow in the entire blade 
channel. 

An investigation of the axisymmetry of the flow was conducted. 

The results are tabulated in Table C-1. The total pressure behind the 
inlet guide vanes, was measured at eight different radial survey 
holes in plane S. P. 2 of Fig. 7. The average reading of 
578.1 counts. The nearest reading to the average was 573.9 counts, 
taken at the hole labeled 33°45*. This hole was used for all readings 
of The relative error was 0.73 per cent. The average of read- 

ings taken behind the third stage stator P^^ was 10, 821.5 counts. 

The nearest reading to the average was 10, 810 counts, taken at the 
hole marked 37°27*30" on Fig. 7. This hole was used for all readings 
of The relative error was 0.10 per cent. Within the measurement 
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capability of the instrument, the flow was considered axisymmetric 
and the most representative readings have been used. 

C-2 Program QNRETA 

The data reduction described in Section 4 was accomplished by 
Program ONRETA (Table C-2). Since the calculations are straight- 
forward and have been described in the body of the thesis, they will 
not be repeated hereo 

The input format is as follows: 



Card 1: 


NCASE 


run number (arbitrary) 


(2I8.4F8.2) 


NRUN 


number of data points per run 




PBAR 


reading of barometer (in. of Hg) 




TEAR 


temperature of barometer (^F) 




RSTAG 


rotor stagger angle (43.8^) 




SSTAG 


stator stagger angle (various 


Card 2; 


TO 


temperature at inlet screen (^F) 


(F8.1,4f8.0 


QO 


velocity head at permanent Pitot- 


F8.1,F8.2) 




static tube (counts) 




PSO 


static pressure with reference to 
atmospheric at permanent Pitot- 
static tube (counts) 




PTIG 


total pressure with reference to atmo- 
spheric ahead of first rotor (counts) 




DELP 


pressure rise across the 3 stages 
(counts) 




TACH 


reading of the speed counter (rpm) 




TRAW 


reading of the torque meter (lb) 



Cards 3 through (NRUN - 1): repeat of Card 2 for remaining data points 



139 



The first five entries on cards 2 through (NRUN - 1) are the arithmetic 
average of four readings taken two minutes apart. The last two entries 
are the average of two readings taken five minutes apart. 

The first section of program output is a printout of the input 
data. The second section is the calculated data for the run. The 
third section is the calculated nondimens ional coefficients. 

Tables C-3 through C-13 present the results of runs 1 through 11. 
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P«riph«ral Angit, 0«gr««« 



TABLE C-1 



INVESTIGATION OF AXISYMMETRY OF FLOW \ 
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Table C-6, Output of Program OKRETA, Run 4, Stator Stagger Angle = 23.8 
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Table C-8. Output of Program ONRETA, Run 6, Stator Stagger Angle = 31.8 
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Table C-9- Output of Program ONRETA, Run 7> Stator Stagger Angle = 31.8 
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Table C-11. Output of Program OKRETA, Run 9> Stator Stagger Angle = 35-8 
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Table C-12. Oubpu't of Program ONRETA, Run 10, Sbabor Sb agger Angle = 39*8^ 
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Table C-13. Output of Program ONRETA, Run 11, Stator Stagger Angle = 44.3^ 



APPENDIX D 



DETAILS OF PREDICTION PROGRAM 

This Appendix contains details concerning AXC03 which may be of 
interest to a limited number of readers. Table D-1 is a listing 
of the main program and the seven subroutines. Section D-1 is a 
description of the method of inserting data. Table D-2 is a sample 
of program output. The rewrite of input data, a sample of interstage 
data, and a sample summary of one case are included. Tables D-3 
through D-16 are the summary printouts of the cases which were compared 
to measured performance of the machine. Table D-17 is a printout of 
program CHECK which verified the errors in the polynomials of sub- 
routines THEORY and CASCAD, 

D-1 Input Data 

There is no limit to the number of cases which may be considered^ 
since data are read in prior to execution for each case, and output 
printed prior to considering the next case. The executive routine 
reads input data and prints it before computation begins. The input 
data consists of thirteen cards per case. 



Card 1: 


NUM 


case number (arbitrary) 


(17, 


GAM 


ratio of specific heats 


4F7.3) 


RN 


referred speed (rpm/V^) 




WREM 


initial referred flow rate 
(lbm"\/^/sec psia) 




DW 


referred flow rate increment 
(lbm“\/^/sec psia) 
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Card 2; 


LMAX 


number of stages (3) 


(17, 


ZETAI 


loss coefficient of inlet duct 


5F7.3) 




(assumed 0.0) 




ZIGV 


loss coefficient of inlet guide vanes 
(assumed 0.05) 




EGV 


efficiency of exit guide vanes (assumed 0.85) 




ETD 


efficiency of diffuser (assumed 1.0) 




EEX 


efficiency of exit (assumed 1.0) 


Card 3: 


DIN 


inlet pipe diameter, in. 


(3F7.3, 


DAV 


average blade diameter, in. 


4f7.1, 


DOU 


outlet pipe diameter, in. 


3F7.3) 


AEA 


2 

annulus area ahead of IGV, in. 




AIGV 


2 

annulus area after IGV, in. 




AGV 


2 

annulus area after EGV, in. 




ADI 


2 

diffuser exit area, in. 




ALE 


inlet guide vane exit angle, degrees 




SMAX 


maximum limit, 29 




SMI 


minimum limit, jq) / ^ j o > ^^8* 29 


Card 4: 


AAA 


an array of annulus areas after rotor one. 


(10F7.3) 




2 

stator one ,.... stator LMAX, in. 


Card 5; 


BH 


an array of blade heights of rotor one...« 


(10F7. 3) 




stator LMAX, in. 


Card 6: 


BLOCK 


an array of blockage factors for rotor one... 


(10F7.3) 




stator LMAX (assumed 1.0) 


Card 7: 


WDO 


an array of work done factors for stage one. . 


(10F7.3) 




stage LMAX (assumed 1.0) 
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Card 8: 


BLDR 


number of blades in the row 




(Rotor 


CAMR 


camber angle, degrees 




data) 


STAGR 


stagger angle, degrees 




(9F7.3) 


SHR 


shape correction factor, (1.0) 






THR 


thickness to chord ratio 






CHR 


blade chord length, in. 






DIAR 


mean blade diameter, in. 






DER 


tip clearance, in 






CDR 


minimum profile drag coefficient 


(assumed 






various) 




Card 9; 


BLDS 


number of blades in the row 




(Stator 


CAMS 


camber angle, degrees 




data) 


STAGS 


stagger angle, degrees 




OF7.3) 


SHS 


shape correction factor (1.0) 






THS 


thickness to chord ratio 






CHS 


blade chord length, in. 






DIAS 


blade mean diameter, in. 






DES 


tip clearance, in. 






CDS 


minimum profile drag coefficienw 


(assumed 






various) 




Cards 10 


and 12 are 


duplicates of 8 




Cards 11 


and 13 are 


duplicates of 9 
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Table D-1. Listing of Program AXC03 (cont.) 
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Table D-1. Listing of Program AXC03 (cont.) 
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Table D-1. Listing of Program AXC03 (cont 
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Table EU2. Sample Output of Program AXC03 



30.000 20.320 42.800 l.OOO O.ICC 2.600 28.800 0.037 0.020 
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CALCULATED AXIAL CCMPRESSOR PERFORMANCE 
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Table X>-2, Sample Output of Program AXC03 (cont.) 



CALCULATED AXIAL COMPRESSOR PERFORMANCE* OVERALL DATA 
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CALCULATED AXIAL COMPRESSOR PERFORMANCE, OVERALL DATA 
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Table D-3. Output of Program AXC03, Case 11, Cj^ = 0.000 
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Table D-4. Output of Program AXC03, Case 12, Crymi n = O.OO 6 
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Table D-5. Output of Program AXC03, Case 13, C- . = 0.008 
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STATOR STAGGER ANGLE = 27 . RO 
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Table D-9* Output of Program AXC03, Case 31 > ^Dmin “ 0.000 
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Table D-17. Program CHECK, Listing and Output (cont.) 



MAIN DATE = 67160 16/30/A8 PAGE 0003 
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Table D-17, Program CHECK, Listing and Output (cont.) 



CCMPARISUN OF AlTLAL VALUtS IN NASA SP-36 WITH VALUES CALCULATED 6Y AXCOM2 



SOLIDITY » 0.8 



LlkO CAMBER INCIDENCE VERSUS ANGLE 

ANGLE, CEGKEES ACTUAL INCIDENCE CALCULATED INCIDENCE ACTUAL - CALCULATED PERCENT ERROR 



0.0 


o.O 


-0.0009329878 


0.0009329878 


0.0932987332 


b. 0000000000 


0.329999<;833 


0.3080989718 


0.0219010115 


6.6366682053 


10.0000000000 


0.6799999475 


0.6169836521 


0.0630162954 


9.2671022415 


lb. 0000000000 


0.9899999499 


0.9257200956 


0.0642798543 


6.4929141998 


/O. 0000000000 


1. 3099994o59 


1.2343072891 


0.0756921768 


5.7780303955 


. b.OOOOOCOUOO 


I,6l999988b6 


1.5427494049 


0.0772504807 


4.76854801 18 


30.0000000000 


1.9299993515 


1.8510446548 


0.0789546967 


4.0909175873 


3b. 0000000000 


2.2299995422 


2. 1591939926 


0.0708055496 


3. 1 751365662 


^O.OuOOOOOOOO 


2.5499992371 


2.4671907425 


0.0828084946 


3.2473926544 


Hb. 0000000000 


2.8299999237 


2.7750444412 


0.0549554825 


1.9418897629 


so.ooouocoooo 


3. 1 199996«56 


3.0827503204 


0.0372495651 


1.1938962936 


b‘). 0000000000 


3.3899993896 


3.3903064728 


-0.0003070831 


-0.0090584978 


to. 0000000000 


3.6199998856 


3.6977157593 


-0.0777158737 


-2.1468467712 


bb.OoOOOOOOOO 


3.8299999237 


4.0049800873 


-0.1749801636 


-4.5686721802 


TO.OOOCOOOOOO 


4.0099992752 


4.3 120946884 


-0.3020954132 


-7.5335502625 



SLOPE FACTOR N VERSUS ANGLE 



t , u tORE E b 


actual factok 


CALCULATED FACTOR 


ACTUAL - CALCULATED 


PERCENT ERROR 


0.0 


-0.0389999971 


-0.0343636572 


-0.0046363398 


11.8880500793 


b.OUOOOOOOOO 


-0.0499999970 


-0.0483931080 


-0.0016068891 


3.2137775421 


lu.ocooooouoo 


-0.0609999970 


-0.0621372499 


0.00 11372529 


-1.8643474579 


ib.OOOOOUOOOO 


-0.0739999413 


-0.0760737062 


0.0020737648 


-2.8023862839 


d’O.OUOOOCOOOO 


-0.0679999995 


-0.0906807184 


0.0026807189 


-3.0462713242 


2t>. JUOOOCOOOO 


-0. 10 39999/23 


-0.1064355373 


0.0024355650 


-2.3418893814 


:: J.OOOOOOOUOO 


-0. 1209999919 


-0. 1238162518 


0.0026162599 


-2.3274869919 


ib.OOUOOOOUOO 


-J. 1409999726 


-0.1433003545 


0.002300381 7 


-1.6314764023 


m;.0000000000 


-0. 1619999409 


-0.1653655767 


0.0033656359 


-2.0775537491 


JOOOOOuOOO 


-0. I8b9999895 


-0. 1904897094 


0.0044897199 


-2.4138278961 


Su.OOuOOOOOOO 


-0.21 19999528 


-0.2191501260 


0.0071501732 


-3.3727235794 


‘)b. OouuOOOOOO 


-0. 2^09999967 


-0.2518250346 


0.0108250380 


-4.4917144775 


c J. OJOOUOGJUO 


-0. 2/69999504 


-0.2889920473 


0.0119920969 


-4.3292751312 


cb.uuUOUOOOOO 


-0. b209999o00 


-0.331 1289549 


0.0101289749 


-3.1554431915 


/C. JOv)OOUOOOO 


-0. 3759999671 


-0.3787130713 


0.0027130842 


-0.7215649486 





OEVIATIUN VERSUS ANGLE 






Cf L.tu^ctS 


AClUAl UcVIATlCN calculated OEVIaTION 


ACTUAL - DEVIATION, 


PERCENT ERROR 


^.0 


o.O 


-0.0014475845 


0.0014475845 


0.1447584033 


0 .OOUOJOOOUO 


0.0599999987 


0.0814747810 


-0.0214747824 


-35.7912902832 


l^.OOJOOOOOUu 


0. 1299999952 


0.1593298316 


-0.0293298364 


-22.5614013672 


1). OOJOOOOOOO 


u. 209999970b 


0.2358062267 


-0.0258062482 


-12.2886896133 


. .OOuOOOJOOO 


0. 2999999523 


0.3145923018 


-0.0145923495 


-4.8641 147614 


^ J . UOUOOOOOOO 


J. 39^9999762 


0.3993784785 


0.0006214976 


0.1553744078 


K-. UOOOOOUCOO 


0. 509999990b 


0.4938524365 


0.0161475539 


3.1661853790 


ib.OOUOUCOOOO 


0.6299999952 


0.601 7042994 


0.0282956958 


4.4913787842 


ft .OOOOUOOOJU 


0. 7599999905 


0.7266232967 


0.0333766937 


4.3916702271 


♦ b . Ut»OCuOOOOO 


0.0999999762 


0.8722926378 


0.0277073383 


3.0785923004 


)O.OOOOliOOCUO 


1 . 0599994659 


1.0424118042 


0.0175876617 


1 .65921^0198 


. 00 JGCMJUUUO 


1 . 2500OOU0UO 


1.2406597137 


0.0093402863 


0.7472229004 


uO.UOOOOOOUOu 


1 .4799995422 


1.4707336426 


0.0092658997 


0.6260744929 


(o .00 JOUOOOO J 


1 . 7500000000 


1.7363204956 


0.0136795044 


0.7816859484 


/O. uCOuuOOuOO 


2.069999t948 


2.0411033630 


0.0288963318 


1.3959579468 
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SLOPE FACTOR M VERSUS ANoLE 



AiNCLLf CEuKEtS 


ACTUAL FACTOR 


CALCULATED FACTOR 


ACTUAL - CALCULATED 


PERCENT ERROR 


0.0 


0.30S9999428 


0.3104162812 


-0.0004163384 


-0.1343027353 


5.0000000000 


0.3119999766 


0.3125761747 


-0.0005761981 


-0.1846788526 


10.0000000000 


0. 3149999976 


0.3147908449 


0.0002091527 


0.0663976669 


15.0000000000 


0.3169999719 


0.3171648979 


-0.0001649261 


-0.0520271473 


20.0000000000 


0. 3209999800 


0.3198031187 


0.0011968613 


0.3728539944 


25.0000000000 


0.3239999413 


0.3228106499 


0.0011892915 


0.3670653105 


io.oooocooooo 


0. 3279999495 


0.3262923956 


0.0017075539 


0.5205957294 


35.0000000000 


0.3329999447 


0.3303530812 


0.0026468635 


0.7948539853 


40.0000000000 


0.3379999995 


0. 3350979686 


0. 0029020309 


0.8585889935 


45.0000000000 


0. 34 39999819 


0.3406320214 


0.0033679605 


0.9790582657 


50.0000000000 


0.3499999642 


0.3470600843 


0.0029398799 


0.8399657607 


55 .0000000000 


0.3369999933 


0.3544869423 


0.0025130510 


0.7039358616 


60.0000000000 


0.3649999499 


0.3630179167 


0.0019820333 


0.5430228710 


^5. 000001.0000 


0.3750000000 


0.3727575541 


0.0022424459 


0.5979855657 


70.0000000000 


0.3879999518 


0.3838108182 


0.0041891336 


1.07967281 34 



ANGLtf LLuKLLS 


SLOPE factor 01/00 VERSUS ANGLE/ 

ACTUAL FACTUK CALCULATED FACTOR ACTUAL - CALCULATED 


PERCENT ERROR 


0.0 


0.0769999623 


0.0754218698 


0.0015780926 


2.0494718552 


b.UOOOUOOOOu 


0.0789999962 


0.0780347586 


0.0009652376 


1.2218189240 


lU.UUOOOOUuOO 


0.0809999704 


0.0802936554 


0.0007063150 


0.8719941378 


15.0000000000 


0.0829999447 


0.0826053023 


0.0003946424 


0.4754729867 


20.COOOOOOOOO 


0.085999965 ( 


0.0853770375 


0.0006229281 


0.7243353128 


25.0JU00C0000 


0.0899999738 


0.0890154839 


0.0009844899 


1.0938777924 


30.0000000000 


0.0939999819 


0.0939274430 


0.0000725389 


0.0771690011 


J5. 0000000000 


0.1019999981 


0.1005203128 


0.0014796853 


1.4506711960 


4U.OCOOOOOOOO 


0.1 10999941b 


0.1092003584 


0.0017995834 


1.6212463379 


45.0000000000 


0. 1209999919 


0. 1203768849 


0.0006231070 


0.5 149644613 


5C. 0000000000 


0.1319999695 


0. 1344547272 


-0.0024547577 


-1.8596649170 


55.00000UJOUU 


0. 1499999 /62 


0. 1518397331 


-0.0018397570 


-1.2265043259 


bO.OOOOOOOOJO 


0. 1719999909 


0.1729413271 


-0.0009413362 


-0.5472884774 


65.0000000000 


0. 1979999542 


0.1981655955 


-0.U001656413 


-0.0836572051 


/O.OQOOOOOOOO 


0.2269999981 


0.2279177904 


-0.0009177923 


-0.4043137431 
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APPENDIX E 



PROPOSED DESIGN OF AN IMPROVED INLET DUCT 



The present inlet duct is deficient because of the shape of its 
bellmouth, because of the small velocities that occur in it during 
operation VThich make accurate determination of volume flow rates 
impossible, and because it has no provision to eliminate effects of 
atmospheric gusts or winds. 

The entrance bellmouth should have more gradual changes in wall 
curvature to eliminate local separations at the entrance to the 
cylindrical duct. It was stated earlier that the velocity head in 
this duct with 36 in. diameter is about 0.6 in. of water, a value too 
small for accurate measurements. To avoid large fluctuations of the 
read-out of the Texas Instruments pressure gauge, for small variation 

in flow rate, the velocity head should be increased to about 6 in. 

2 3 

of water or 33 Ib/ft , at a volume flow rate of 360 ft /sec, giving 

a response at the instrument of about 8000 counts. For an assumed 

2 

incompressible flow with a mass density ^ ^ of 0.002378 lb sec / 

4 2 

ft , and a velocity head of 33 Ib/ft , the velocity in the 

measuring plane of the inlet duct would have to be 









(ft/sec) 



( 6 ) 



2 

Hence the required flow area is 360/166 =2.16 ft , corresponding 
to a circular area with 19.9 in. diameter. Because of the small 
pressure rise of the compressor it is necessary to convert the 
velocity head at the measuring plane into static pressure rise by 
arranging an optimum diffuser between this section and the inlet 
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pipe to the compressor that has a diameter of 36 in. Test data of 
Reference 11 show that the included diffuser angle should be 8^ for 
the present application. Larger diffuser angles produce separation 
with associated flow instabilities, and smaller angles require ducts 
with increased lengths where the boundary layer growth is excessive 
and reduces pressure recovery. For a diffuser angle of 8^ the length 
L of the diffuser is 



L 



- n.‘) 

Z -fcan 4 ° 



\ I G f n* 



Figure E-1 shows the design of the proposed inlet duct. A well 
designed bellmouth from a diameter of 48 in. to the throat of 19.9 
in. requires an axial length of about 26 in. , giving a total length 
of the duct of 116 + 26 = 142 in. Ahead of the bellmouth inlet a 
structure will be attached to support several layers of fine mesh 
screen to eliminate flow disturbances by gusts in the surrounding 
atmosphere. The whole duct is mounted on a trolley which can be 
rolled onto the apron outside the test cell to permit simple attach- 
ment to and removal from the presently installed cylindrical inlet 
duct. It is recommended also that an additional honeycomb flow 
straightener be installed at the entrance of this duct. 

The bellmouth and the diffuser could be molded out of plastics, 
reinforced by fiberglass, by using a wooden template which is split 
at the smallest diameter. 
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